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ABSTRACT 


Title:  Irreversible  Thermodynamic  Effects  in  Inhomogeneous 

Media  and  Their  Applications  in  Certain  Geoelectric 
Problems . 


Author:  Bijan  Nourbehecht 


Submitted  to  the  Department  of  Geology  and  Geophysics 
on  November  19,  1962  in  partial  fulfillment  of  the 
requirements  for  the  degree  of  Doctor  of  Philosophy 
at  the  Massachusetts  Institute  of  Technology 


A  theoretical  Investigation  of  irreversible  coupling 
effects  in  the  earth  is  carried  out  to  get  a  better  under¬ 
standing  of  the  role  of  coupling  in  producing  spontaneous 
polarization  currents.  It  is  found  that  electrokinetic 
and  electrochemical  couplings  are  the  most  important,  and 
in  certain  cases,  the  thermoelectric  coupling  can  also  be 
% important.  The  electric  current  is  found  to  depend  on  the 
gradient  of  a  potential  -f  >  called  the  "total  electric 
potential,"  which  satisfies  the  Laplace's  equation.  The 
boundary  value  of  j-  depends  on  the  driving  pressure  P  , 

(or  chemical  potential  J//.  )  and  the  electric  potential p  . 

The  effect  of  geometry  is  investigated  by  solving  the 
equation  for  using  horizontally  and  vertically  layered 
earths,  buried  slab,  and  half  space  with  random  coupling 
properties  as  models.  It  is  found  that  when  the  equal 
pressure  surfaces  cut  across  the  rock  boundaries,  an 
electric  current  is  generated  as  the  surface.  The  magni¬ 
tude  of  this  current  depends  on  the  geometry  as  well  as 
the  differences  in  coupling  properties. 

The  ground  potentials  observed  following  an  under¬ 
ground  nuclear  detonation  are  shown  to  be  induced  by  electro- 
kind  tic  coupling.  The  computed  streaming  potentials  for 
one  layer  over  a  half  space  are  found  to  agree  fairly  well 
with  the  field  data.  Quantitative  interpretation  of  the 
results  gave  a  depth  estimate  slightly  higher  than  the 
actual  depth. 

Electrochemical  diffusion  is  found  to  be  capable  of 
producing  only  small  S.  P.  anomalies  of  a  few  tens  of 
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millivolts,  because  the  variations  in  the  concentrations  of 
the  diffusing  ions  are  small.  However,  the  diffusion 
potentials  are  of  sufficient  magnitude  to  account  for  much 
of  the  background  noise.  In  the  presence  of  electronic 
conductors,  the  self-potentials  are  determined  by  the  con¬ 
centrations  of  reacting  species  (e.g.  FZf+  and  7T ’ttt)  .  When 
iron  is  the  only  important  reacting  species,  the  maximum 
potentials  are  of  the  order  of  one  volt.  The  effect  of 
geometry  on  the  rms  value  of  the  background  noise  is  evalu¬ 
ated  by  using  a  model  composed  of  randomly  distributed 
vertical  dikes. 
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CHAPTER  I 


INTRODUCTION 


Purpose 

For  some  thirty  years ,  thermodynamics  of  irrevers¬ 
ible  processes,  as  developed  by  Onsager  and  others,  has 
been  applied  successfully  to  obtain  important  results  for 
systems  that  are  not' in  equilibrium  and  where  there  is 
an  interaction  of  several  processes.  When  two  or  more 
flows  (e.g.  heat  flow,  electrical  current,  etc.)  are 
taking  place  simultaneously  in  the  same  system  under  the 
action  of  thermodynamic  forces,  such  as  temperature  or 
concentration  gradients,  the  rate  of  each  flow  is  regarded 
as  being  proportional  not  only  to  the  corresponding 
force,  but  also  it  is  supposed  to  be  dependent  on  the 
other  forces  as  well.  These  phenomena,  known  as  "coupling" 
of  flows,  give  rise  to  several  interesting  effects. 

Among  them,  thermoelectricity,  streaming  potential,  and 
electro-osmosis  are  well-known. 

Here  we  shall  apply  the  thermodynamics  of  irrevers¬ 
ible  processes  to  the  study  of  coupling  of  electrical 
current  with  solvent  flow,  heat  flow,  or  ionic  flow  which 
are  caused  by  pressure,  thermal,  or  concentration  gradients. 
In  particular,  we  shall  apply  the  theory  to  a  detailed 
study  of  the  electrical  after-effects  of  an  underground 
nuclear  blast  and  the  currents  generated  across  a  dike 
or  the  contact  between  two  media  in  the  presence  of  a 
chemical  potential  gradient.  This  approach  has  the  advan¬ 
tage  of  unifying  several  apparently  unrelated  phenomena 
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and  giving  quantitative  results  in  each  case  once  the 
relevent  parameters  of  the  medium  are  measured  experi¬ 
mentally. 

Underground  Nuclear  Blast 

The  first  objective  of  this  thesis  is  a  theoretical 
study  of  the  electrical  aftereffects  of  a  buried  nuclear 
explosion  and  the  manifestation  of  these  effects  at  the 
surface.  It  is  possible  that  measurements  of  blast¬ 
generated  surface  potentials  could  be  used  for  an  on-site 
inspection  procedure. 

One  of  the  important  aftereffects  of  an  underground 
blast  is  the  release  of  a  large  amount  of  heat.  A  one 

TO 

kiloton  bomb  releases  about  5  X  10  joules  (Johnson 
and  al.,  1959) ,  and  assuming  a  specific  heat  of 
.2  cal/gm-°K  and  a  density  of  rocks  of  2.5  gm/cm^,  it 
is  easy  to  see  that  a  20KT  bomb  can  heat  to  the  boiling 
point  temperature  of  water  (100°Q)  all  the  rocks  in  a 
sphere  of  175  feet  radius  and  originally  at  20°C.  But 
if  the  explosion  is  totally  contained,  even  though  the 
released  energy  is  enormous,  no  appreciable  temperature 
rise  will  be  observed  at  the  surface  in  a  reasonable 
time  of  months  or  even  years  due  to  the  poor  conductivity 
of  the  rocks.  The  dimensionless  parameter  (Carslaw  and 
Jaeger,  19^7)  for  heat  conduction  is 


where  K  is  the  thermal  diffusivity  and  for  rocks  has  a 
value  of  about  0.01  cm2 /sec;  Y*  is  distance  in  cm;  and 
t  is  time  in  seconds.  For  distances  of  the  order  of 
100  meters  or  more,  and  expressing  time  in  days,  the 
heat  diffusion  parameter  6  is 

©  ^  lO'^t 


(t  in  days) 
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and  therefore  at  these  distances,  a  time  scale  of  years 
is  involved.  Consequently  for  all  practical  purposes, 
the  thermal  disturbances  are  completely  masked  at  the 
surface  (for  an  entirely  contained  explosion)  and  can 
not  be  used  directly  to  detect  the  presence  of  a  blast. 
However,  studies  by  Madden  and  others  (Madden  and  Marshall, 
1959;  Madden  and  Cantwell,  1962;  see  also  table  2.4) 
have  shown  that  rocks  can  couple  efficiently  fluid, 
thermal,  and  chemical  potential  gradients  into  electric 
potentials.  Thermoelectric  coupling  coefficients  for 
rocks  can  be  as  high  as  1  mv/°C  with  average  values  of 
about  .2-. 4  mv^C.  Therefore,  a  temperature  difference 
of  50  to  100°C  can  give  rise  to  a  very  appreciable 
potential  difference. 

Another  important  aftereffect  connected  with  the 
release  of  energy  is  a  pressure  build-up.  Some  of  the 
released  energy  is  used  up  mechanically  to  fracture  the 
rocks  or  produce  seismic  waves;  but  the  greater  part 
of  the  energy  goes  into  heating  up  the  rocks  and  turn¬ 
ing  their  pore  fluids  into  steam.  How  large  the  pressure 
builds  up  and  how  long  it  lasts,  depends  on  several 
factors  such  as  the  size  of  the  bomb,  the  fluid  content 
of  the  rocks,  the  permeability,  the  heat  dissipation, 
and  the  extent  of  fracturing.  If  the  blast  is  wholly 
contained  and  the  fracturing  does  not  open  up  large 
cracks,  the  only  way  this  excess  pressure  can  dissipate 
is  through  diffusion  of  heat  and  condensation  of  the 
stream,  or  outward  flow  of  fluids.  The  first  of  these 
two  mechanisms  is  probably  unimportant,  due  to  the  slow 
rate  of  heat  diffusion  in  rocks.  But  outward  flow  of 
fluids  is  expected  to  be  a  very  important  cause  of 
pressure  decrease  in  many  shots.  The  rate  of  flow 
depends  on  the  permeability  coefficient  which  can  vary 
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wildly  from  one  situation  to  another  by  as  much  as  seven 
orders  of  magnitude  (see  Table  2.3).  Actually,  the 
fracturing  opens  up  many  cracks,  and  increases  the 
effective  permeability.  It  is  possible  that  in  certain 
cases,  cracking  of  the  rocks  opens  up  fissures  reaching 
the  surface  and  thus  allowing  all  of  the  steam  to  escape 
into  the  air.  Even  in  these  cases,  we  can  expect  a 
pressure  effect,  because  the  pressure  inside  the  cavity 
will  be  equal  to  the  atmospheric  pressure  while  the 
ambient  hydrostatic  pressure  will  be  much  higher.  A 
backflow  will  then  take  place.  For  a  sho.tlburied  at  a 
depth  of  200  meters,  this  backflow  pressure  is  of  the 
order  of  tens  of  atmospheres. 

Due  to  all  the  various  factors  mentioned,  it  is 
extremely  difficult  to  predict  generally,  the  behaviour 
of  pressure  with  time  for  any  particular  situation,  but 
we  can  nevertheless  get  a  rough  idea  of  the  time  scales 
involved  by  using  the  following  model.  Assume  that  all 
the  steam  Inside  a  spherical  cavity  of  radius  Ya  has  been 
released  to  the  air  through  fractures  that  reach  the 
surface,  and  fluid  starts  flowing  into  the  cavity  under 
the  action  of  a  uniform  hydrostatic  pressure  P„  (we  are 
neglecting  the  asymmetry  due  to  gravity).  This  model, 
although  very  simple,  will  give  an  idea  of  the  time 
scales  involved,  and  is  expected  to  be  useful  in  certain 
cases . 

The  flow  of  incompressible,  fluid  inside  the.  earth 
is  governed  by 


O'  =  -  K  VP 


1.2 


VZF  =  0 


1.3 
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where  J  is  the  rate  of  flow  and  K  is  the  permeability 
of  the  medium.  Solution  to  1.3  having  spherical 
symmetry  is 


F  -  P(i  -  £)  1.4 

where  is  the  hydrostatic  pressure  at  large  distances 

from  the  cavity.  The  total  flux  across  the  cavity 
boundary  is 

Q.  =  If  (rMA)  =  -4TTKr;2(^)r#  1.5 

Substituting  from  1.4  into  1.5,  we  get  for  the  rate  of 
flow  of  fluid  inside  the  cavity 

,  1 

1.6 

and  the  time  it  takes  to  fill  up  is  approximatively 

taV_  _  .  I*  j  , 

q  4tiKr,F, 

Taking  some  typical  values  of  =  50  meters,  10 

atmospheres,  and  K  «■  10~"^MKS,  we  get 

1  zsoo  n 

t  =  - — - -  £  Z.\  X\t>  Me.  =  Z^O  dcuiS 

This  time  constant  can  vary  easily  by  a  factor  of  100, 
and  the  estimate  given  here  is  probably  too  high  because 
we  have  neglected  the  effect  of  fracturing.  This 
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calculation  was  done  to  show  that  time  scales  of  a  few 
days  to  even  years  are  possible  for  the  electrokinetic 
effect . 

Rocks  have  rather  large  electrokinetic  coefficients 
of  the  order  of  a  few  millivolts  per  atmosphere — some¬ 
times  as  high  as  30  mv/atm — and  therefore  a  pressure 
difference  of  only  ten  atmospheres  is  sufficient  to 
give  potential  differences  of  50  millivolts  or  better. 

Other  after-effects  such  as  gravity  anomaly  or  seismic 
disturbances  are  not  discussed  here  as  they  have  no 
direct  bearing  on  this  thesis.  A  list  of  all  posible 
effects  and  a  brief  discussion  of  the  importance  of 
each  effect  is  given  by  T.  R.  Madden  (Madden  and  Cantwell, 
I960;  p.7).  In  many  ways  the  nuclear  detection  problem 
is  very  similar  to  mining  geophysics  problems  because  of 
the  geometry  and  size  of  the  target,  and  therefore  con¬ 
ventional  methods  based  on  classical  effects  (gravity, 
resistivity,  etc.)  could  be  used.  But  it  differs  from 
mining  problems  because  of  the  importance  of  streaming 
potential  and  thermoelectric  phenomena.  Mineralized 
areas,  due  to  the  rather  long  time  scales  Involved,  have 
reached  equilibrium  at  least  as  far  as  temperature  and 
pressure  are  concerned;  but  a  blast  area  is  far  from 
having  reached  equilibrium  and  therefore  these  pressure 
and  temperature  effects  play  an  important  role  in  pro¬ 
ducing  natural  currents. 

Spontaneous  Polarization 

The  second  objective  of  this  thesis  is  a  theore¬ 
tical  and  quantitative  study  of  the  causes  of  certain 
types  of  spontaneous  polarization  currents.  Spontan¬ 
eous  polarization  effect  or  S.  P. ,  which  can  be  produced 
by  many  different  causes,  are  well-known  and  many  authors 
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have  discussed  qualitatively  several  of  the  mechanisms 
responsible  for  these  potentials  (see  Dobrin,  1952, 
p.289  and  378;  Heiland,  1946,  p.668;  Mounce  and  Rust, 

194-5;  Wyllie,  1949)’  Most  of  these  explanations  depend 
on  oxidation-reduction  reactions  or  concentration  cells. 
One  of  the  more  quantitative  explanations  has  been 
advanced  by  Sato  who  shows  that  the  S.  P.  currents 
around  a  sulfide  body  are  intimately  connected  with  ore 
body,  and  that  spontaneous  polarization  could  take  place 
even  if  the  ore  body  itself  did  not  undergo  any  changes 
and  acted  only  as  an  inert  conductor  (Sato  and  Mooney, 
i960).  Others,  studying  the  causes  of  S .  P.  in  drill 
holes  (Wyllie,  1949),  attribute  the  polarization  to  the 
presence  of  a  shale  layer  which  acts  as  an  almost 
perfect  semi-permeable  membrane  and  stops  all  negative 
ions  from  diffusing  through,  while  it  allows  the  positive 
ions  to  diffuse  almost  freely. 

The  methods  of  irreversible  thermodynamics  are  used 
here  to  unify  all  the  various  explanations  based  on 
electrochemical  or  electrokinetic  couplings.  Shale 
barriers  represent  an  extreme  case  and  may  give  the 
largest  potentials,  but  are  not  lndispensible  for 
causing  S.  P.  currents.  Any  two  materials  in  contact, 
with  dissimilar  diffusion  rates,  will  tend  to  produce 
S.  P.  if  there  is  a  chemical  potential  gradient  cutting 
across  their  common  boundary.  Oxidation-reduction 
problems  are  dealt  with  in  much  the  same  way  as  concen¬ 
tration  diffusion  problems.  In  particular  the  effects 
of  geometry  can  be  handled  very  nicely  using  this  approach 
Specifically,'  we  have  investigated  the  effect  of  vertical 
contacts  and  dikes  on  the  magnitude  of  the  S.  P.  anomalies 

Historical  Review 


The  original  ideas  of  this/  thesis,  which  are  based 
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on  principle  of  irreversible  thermodynamics,  were 
recently  suggested  by  T.  R.  Madden  in  connection  with 
the  possibility  of  using  electric  methods  for  onsite 
detection  of  underground  nuclear  blasts  (Madden  and 
Cantwell,  i960.)  Because  of  the  non-equilibrium  con¬ 
ditions  in  the  near  vicinity  of  the  blast,  and  the 
large  coupling  coefficients  of  the  rocks,  it  was  thought 
that  measurements  of  electrical  aftereffects  caused  by 
the  coupling  of  fluid  and  thermal  gradients  with 
electric  potentials  could  be  used  to  detect  large 
underground  nuclear  blasts.  However,  it  can  be  assumed 
that  these  gradients  are  confined  to  depth  and  vanish  on 
the  surface,  and  it  is  therefore  necessary  to  undertake 
a  theoretical  study  of  the.  problem  to  find  out  the 
surface  expression  of  the  effects  which  are  produced  in 
the  near  vicinity  of  the  source.  As  the  research  got 
underway,  it  became  obvious  that  these  same  principles 
could  be  used  to  explain  some  of  the  causes  of  spontan¬ 
eous  polarization  currents,  and  the  effect  of  geometry 
on  their  magnitudes. 

Many  of  these  coupling  phenomena  have  been  known 
for  some  time  but  their  description  based  on  thermodynamics 
of  irreversible  processes  ds  rather  recent.  Onsager 
developed  his  theory  of  irreversible  thermodynamics  in 
1931  and  applied  it  to  the  study  of  chemical  monomole- 
cular  triangle  reactions  and  heat  flow  in  anisotropic 
crystals  (DeGroot,  1961;  Denbigh,  1951) •  Later  on,  the 
theory  was  applied  by  DeGroot  and  others  to  describe  all 
kinds  of  coupling  of  flows  such  as  electrochemical  and 
thermoelectric  effects. 

Schlumberger  and  Leonardon  (193^a  and  193^6)  studied 
the  causes  of  S.  P.  potentials  in  bore  holes  and 
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concluded  that  they  are  due  to  electrofiltration  and 
electrochemical  couplings.  They  even  made  some  labora¬ 
tory  studies  that  supported  their  view.  But  their 
theoretical  explanation  was  rather  qualitative  and  based 
on  Nerst's  equations  for  e.m.f.  of  reversible  cells. 

These  equations  are  of  limited  application  in  most  cases 
of  spontaneous  polarization.  Most  other  authors  who 
have  studied  this  problem,  have  systematically  followed 
Schlumberger  in  explaining  the  cause  of  S.  P.  in  bore¬ 
holes  . 

Mounce  and  Rust  (19^5)  studied  the  same  problem 
and  showed  experimentally  that  currents  can  be  set  up 
in  a  system  composed  of  two  aqueous  salt  solutions  of 
different  concentrations  and  a  shale  layer  which  acts 
as  a  semi-permeable  membrane ;  They  supposed  that  a 
shale  layer  is  essential  to  the  creation  of  currents. 

Wyllie  (19^9)  took  up  the  same  problem  but  he  came 
up  with  no  new  results.  He  discussed  the  role  of  the 
shale  layer  and  stated  that  "action  of  a  shale  barrier 
is  analogous  to  a  glass  membrane  separating  two  acid 
solutions  of  different  hydrogen  ion  concentrations. 

The  shale  behaves  as  a  sodium  electrode  and  is  responsive 
to  the  activities  of  the  sodium  ions  in  the  two  solutions 
in  such  a  way  that  the  potential  can  be  calculated  by 
means  of  the  Nerst  equation."  The  shale  membrane  is  an 
extreme  case  because  it  stops  ions  of  one  sign  and  only 
lets  through  ions  of  the  other  sign.  Actually  for 
creating  electrochemical  potential  it  is  only  necessary 
to  have  an  asymmetry  in  the  diffusion  rates.  Wyllie 
(1955)  also  underestimated  the  effect  of  geometry  and 
inhomogeneities  on  S.  P.  potentials  in  drill  holes  and 
writes  "These  effects,  though  serious  in  electrical  log 
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interpretation,  are  purely  physical.  They  may  be  largely 
overcome. . .by  the  application  of  static  (bucked  current) 

S.  P.  logging  techniques..." 

Cagniard  (1956)  considered  the  role  of  electro¬ 
filtration  as  a  source  of  spontaneous  polarization.  He 
states  that  this  phenomenon  is  the  cause  of  observed 
negative  anomalies  associated  with  topographic  highs. 

Kermabon  (195 6)  made  some  measurements  of  the 
electrokinetic  coupling  and  permeability  coefficients 
of  rocks  and  found  that,  in  general,  rocks  have  strong 
coupling  coefficients. 

Madden  and  Marshall  (1959)  made  a  systematic  study 
of  electrokinetic,  thermoelectric,  and  electrochemical 
phenomena  in  connection  with  causes  and  background 
effects,  of  Induced  polarization  phenomenon.  To  the 
author's  knowledge  this  is  the  first  application  of 
Onsager's  theory  to  geophysical  problems. 

Since  the  importance  of  coupling  phenomena  in  many 
geophysical  problems  has  become  well  established.  Madden 
and  Cantwell  (1962)  have  undertaken  an  extensive  program 
of  laboratory  measurements  of  electrokinetic,  thermo¬ 
electric,  and  electrochemical  coefficients  of  rocks. 

Outline  of  the  Thesis 

The  Laboratory  measurements  show  how  efficiently 
rocks  can  couple  fluid,  thermal,  and  chemical  gradients 
into  electric  potentials.  This  investigation  was  under¬ 
taken  (;o  understand  better  how  these  coupling  phenomena 
which  are  confined  to  the  near  vicinity  of  the  sources — 
where  the  gradients  exist — manifest  themselves  on  the 
surface,  and  what  effects  rock  inhomogeneities  have  on  the 
magnitudes  of  these  couplings. 
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In  the  first  part  of  Chapter  I,  the  role  of 
electrokinetic  and  thermoelectric  couplings,  produced  by 
an  underground  nuclear  blast,  are  discussed  in  view  of 
their  importance  in  producing  electric  currents  around 
the  blast  area.  It  is  shown  that  the  time  scales  involved 
are  always  of  the  order  of  years  for  the  thermoelectric 
effect  and  can  be  of  the  order  of  years  for  the  electro- 
kinetic  effect.  Furthermore,  electrochemical  and  electro- 
kinetic  couplings  are  discussed  as  possible  causes  of 
the  spontaneous  polarization  currents  associated  with 
sulfide  ore  bodies . 

In  Chapter  II,  Onsager's  phenomenological  relations 
are  introduced  without  any  attempt  at  justifying  them. 

The  connection  between  certain  of  the  phenomenological 
coefficients  and  more  familiar  electrochemical  and  thermal 
quantities  are  discussed.  Starting  from  Onsager's 
relations,  it  is  shown  that  the  same  set  of  equations  and 
boundary  conditions  govern  the  electrokinetic,  thermo¬ 
electric,  and  electrochemical  couplings  in  the  earth.  The 
effect: of  the  driving  pressure  (or  temperature,  chemical 
potential,  etc.)  is  shown  to  be  equivalent  to  that  of 
electric  dipole  sources  distributed  along  the  boundaries 
of  the  different  rock  zones . 

In  Chapter  III,  the  equation  governing  the  coupling 
of  flows  is  solved  for  some  simple  geometries.  The 
geometry  is  shown  to  have  an  essential  role  in  creating 
spontaneous  polarization  currents  detectable  at  the 
surface.  It  Is  shown  that  surface  electrical  signals 
arise  whenever  inhomogeneities  in  rock  properties  cut 
across  the  equl-pressure  (or  temperature  or  chemical 
potential)  surfaces  associated  with  the  flow  of  fluid  (or 
heat  or  ions).  The  relevent  parameter,  on  which  the 
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magnitude  of  the  electric  potential  depends,  is  the 
ratio  of  the  coupling  coefficient  of  the  rock  to  its 
electrical  conductivity.  Several  geometries  of  horizon¬ 
tally  and  vertically  layered  earths  are  used  to  compute 
the  thermoelectric,  the  electrokinetic,  and  the  electro¬ 
chemical  potentials  observed  on  the  surface.  The  effect 
of  geometry  on  the  background  S.  P.  noise  is  computed 
for  a  particularly  simple  model  of  vertically  distri¬ 
buted  random  variations  in  rock  properties. 

Chapter  IV  is  devoted  to  the  applications  of  these 
results  to  the  bomb  detection  problem.  Maximum  potentials 
to  be  expected  from  thermoelectric  and  electrokinetic 
couplings  are  computed  to  get  the  upper  bound  for  these 
effects.  It  is  shown  that  electrokinetic  coupling  is 
probably  the  most  important  cause  of  the  potentials 
associated  with  underground  nuclear  blasts,  while 
thermoelectric  coupling  might  reinforce  these  potentials 
in  certain  cases.  Computed  magnitudes  of  the  potentials 
are  of  the  same  order  of  magnitude  as  potentials 
measured  during  the  Nevada  underground  test  shots. 

Chapter  V  is  concerned  with  the  spontaneous  polari¬ 
zation  currents  associated  with  sulfide  ore  bodies, 
dikes,  and  other  inhomogeneities.  Electrochemical 
coupling  is  believed  to  be  the  most  important  cause  of 
these  currents  at  least  for  small  scale  effects,  with 
thermoelectric  and  electrokinetic  couplings  may  be 
important  in  certain  cases.  Coupling  phenomena  are 
examined  for  their  role  in  producing  S.  P.  noise. 

Chapter  VI  is  devoted  to  a  brief  discussion  of  some 
of  the  problems  in  electrochemistry  and  chemical  properties 
of  pore  fluids  which  emerge  from  the  investigations  of 
this  thesis. 


CHAPTER  II 


REVIEW  OP  IRREVERSIBLE  THERMODYNAMICS 


In  Chapter  I,  it  was  stated  that  the  thermodynamics 
of  irreversible  processes  can  be  used  advantageously  to 
study  a  whole  host  of  apparently  unrelated  phenomena. 
Spontaneous  polarization  and  electric  currents  generated 
by  thermal  gradients  in  non-mineralized  areas  are  examples. 

In  section  I  of  this  chapter,  Onsager's  theory  of  irrever¬ 
sible  thermodynamics  is  introduced  in  a  very  qualitative 
way  with  emphasis  on  the  meaning  of  the  theory  rather  than 
the  origin  and  proof  of  the  theory. 

The  remaining  sections  of  this  chapter  are  concerned 
with  the  forces  and  flows  in  Onsager's  theory  and  their 
application  to  study  of  coupling  effects  in  the  earth. 

Section  I  -  Onsager's  Relations 

Classical  thermodynamics — which  could  more  properly 
be  called  thermo-statics — deals  with  states  of  equilibrium, 
and  modern  irreversible  thermodynamics  as  developed  by 
Onsager,  DeGroot,  and  others  deals  with  states  of  non¬ 
equilibrium,  when  the  system  under  consideration  is  dis¬ 
placed  from  equilibrium  and  flows  are  taking  place.  This 
theory  has  had  its  greatest  successes  in  the  study  of 
steady-state  time  Invariant  systems,  but  it  has  been  also 
useful  in  the  study  of  slowly  varying  systems  for  small 
flow  densities. 

By  analogy  with  dynamics,  Onsager  regards  the  processes 
of  diffusion,  heat  flow,  matter  flow,  and  current  flow  as 
being  representable  in  terms  of  velocities  and  "thermody¬ 
namic  forces"  which  are  the  negative  gradients  of  temperature 
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of  chemical  potential,  etc.,  and  depend  on  how  far  the 
system  is  displaced  from  equilibrium.  For  example,  in 
the  flow  of  electricity,  the  velocity  is  equal  to  the 
current  density  J  and  the  force  is  the  negative  gradient 
of  the  electric  potential  -v^  .  When  these  flows  take 
place  separately,  they  are  governed  by  the  well-known 
empirical  relations 


3^=  -  KVT  (Fourier's  Law)  2.1 

3S  __  —  K'VP  (Darey's  Law)  2.2 

3^  =  —liVC  (Fick's  Law)  2.3 

3  =■  -<5 ■'7^  (Ohm's  Law)  2.4 


where  n  „  T  3  and  3  are  flows  of  heat,  solvent, 
matter,  and  current;  VT;  VP,  VC  and  V0  are  gradients 
of  temperature,  pressure,  concentration,  and  electric 
potential;  K.  )  D  and  CT  are  generalized  conductiv¬ 
ities.  Therefore,  in  general,  if  3  is  the  rate  of 
flow,  X  is  the  force,  and  L  is  a  scalar  quantity,  there 
is  an  empirical  linear  relationship  of  the  form 


3  =  LX  ^  2-5 

where  Lj  has  the  significance  of  "conductivity." 

When  two,  three,  or  more  flows  take  place  simultan¬ 
eously,  there  is  interference  and  coupling  between  these 
flows.  Peltier  effect  (evolution  of  heat  at  junctions 
of  metals  due  to  flow  of  an  electric  current)  and  the 
thermoelectric  force  (e.m.f.  resulting  from  the  mainten¬ 
ance  of  the  junctions  at  different  temperature)  are 
well-known  examples  of  such  coupling  between  two  flows. 
When  several  processes  take  place  simultaneously. 
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equations  of  the  form  2.5  are  inadequate  to  describe  the 
flows  and  In  the  absence  of  proof  to  the  contrary,  it 
must  be  assumed  that  each  flow  depends  on  all  the  forces 

3;  -  5  (x., *>,•••'  Xh)  2-6 


For  a  complete  discussion  of  these  coupling  effects,  the 
reader  is  referred  to  Callen  (i960),  DeGroot  (1961), 

Denbigh  •  (1951) *  Fitts  (1962)  and  Prigogine  (1955)*  Only 
the  most  important  results  are  listed  here. 

The  first  postulate  of  thermodynamics  of  irreversible 
processes  is  that  in  a  system  in  which  several  flows 
are  taking  place  simultaneously  due  to  the  action  of  several 
forces  X-  >  it  is  possible  to  express  each  flow  as  a 
linear  combination  of  all  these  forces: 

*) 

These  equations  are  called  the  phenomenological  relations 
or  thermodynamic  equations  of  motion.  Remembering  that 
in  the  absence  of  all  forces,  the  flows  are  zero,  it  is 
easy  to  see  that  the  equation  2.7  is  a  Taylor's  series 
expansion  of  equation  2.6  approximated  to  the  first  linear 
term.  The  assumption  of  linearity  is  by  no  means  obvious, 
and  needs  to  be  verified  in  each  case.  In  general,  it 
is  valid  when  the  system  is  displaced  only  slightly  from 
equilibrium,  but  terms  of ' higher  orders  should  be  included 
when  the  system  is  far  from  equilibrium.  An  approach  to 
include  non-linear  terms  is  indicated  in  Prigogine  (1955* 
p.93).  It  seems,  however,  that  except  in  the  study  of 
rates  of  chemical  reactions,  one  is  never  too  far  from 
equilibrium  and  therefore  the  linear  approximation  is  a 
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good  one.  Equation  2.7  states  that  each  flow  is  caused 
by  contribution  of  all  forces  through  the  phenomenoligical 
coefficients  by  .  The  coefficients  are  for  example 
the  heat  conductivity,  the  ordinary  diffusion  coefficient 
and  the  electrical  conductivity.  The  coefficients  L»4j' 
with  C  ^ j  are  the  coupling  or  "drag"  coefficients  (DeGroot, 
1961). 

The  second  postulate  is  the  principle  of  microscopic 
reversibility  which  is  at  the  basis  of  Onsager's  theorem. 
For  a  complete  discussion  of  the  principle  of  miscroscopic 
reversibility,  .the  reader  is  referred  to  DeGroot  (1961), 
Denbigh  (1951)  and  Prigigine  (1955).  Onsager's  theorem 
states  that:  If  the  forces  and  flows  are  chosen  in  such 
a  way  that  the  sum  of  the  products  adds  up  to  the 

rate  of  creation  of  entropy  multiplied  by  the  temperature, 
the  matrix  of  phenomenoligical  coefficients  is 
symmetric 

If  Tg-  £.8 

Then  L»,y  3  by*,'  2.9 

This  result  implies  that  there  is  a  symmetry  in  the 
coupling  of  the  various  processes.  The  choice  for  these 
flows  and  forces  will  be  indicated  later  on.  Onsager's 
relations  Ly»  by- which  have  been  proved  in  the  case  of 
fluctuations  about  equilibrium,  are  assumed  to  be  valid 
in  case  of  irreversible  flows  for  small  flow  densities. 
Onsager's  theory  does  not  indicate  whether  or  not  the 
thermodynamic  equations  of  motion  are  correct.  In  each 
case  it  has  to  be  found  by  experiment  whether  or  not 
these  equations  are  a  good  approximation,  and  in  parti¬ 
cular  it  has  to  be  shown  that  the  forces  V.  which 
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satisfy  equations  2.7  and  2.8  are  those  to  which  the  flows 
are  in  practice  proportional. 


Section  II  -  Forces  and  Flows  in  Onsager's  Relations 

The  two  postulates  of  irreversible  thermodynamics, 
leave  a  great  latitude  for  the  choice  of  flows  and  forces; 
as  a  matter  of  fact,  there  are  differences  among  various 
authors  in  the  choice  of  the  various  flows  and  forces. 

Here  we  shall  follow  Denbigh  as  modified  by  Eckardt 
and  Callen  (DeGroot,  1961,  p.28;  Madden  and  Marshall, 

1951,  p.23).  For  the  processes  which  are  of  interest  to 
us: 


Heat 

Solvent 

Matter 


4 


Flow 

Jq 


Force 


p 

calories/cm  -sec 

cm^/cm2-sec 

modes/cm2-sec 


-vr 

-VP 

-  VJU-  FiV<p  ( 

(non-hm’c) 


=  electrical  potential 
T  =  temperature 
P  =  pressure 


J~i  =  chemical  potential 
F  =  Faraday's  constant 
2  =  ionic  valence 


The  current  flow  is  not  shown  explicitly  in  the  above 
equations,  because  in  electrolytic  solutions,  the  current 
flow  is  not  independent  of  the  flow  of  matter,  but  is  a 
lirlear  combination  of  the  anionic  and  cathionic  flows 


1  -  F  H  **  i  •  •) 


2.10 


In  the  above  ..equations,  we  have  followed  Eckardt  in 
defining  the  "heat  flow"  vector  (DeGroot,  1961)  as 
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where  3E  is  the  total  "energy  flow"  and  is  the  matter 
flow.  In  the  absence  of  matter  flow  the  distinction  between 
heat  flow  and  energy  flow  is  irrelevent. 

Assuming  that  the  matter  flow  consists  of  only  one 
cathionic  flow  and  one  anionic  flow  ,  we  can  write 

Onsager's  relations  in  the  following  matrix  form 


r  i 

^*i»  4.^ 

-7T/T 

3a 

U}1  4*,^  Ljj.^ 

-  VP 

4»}»  i-jj 

3« 

.Uii 

The  extension  to  the  more  general  case  of  many  ionic  flows 
is  straight-forward.  From  equations  2.11,  the  possible 
magnitudes  of  coupling  effects  can  be  deduced  in  terms  of 
the  coupling  coefficients L>,j  .  As  an  example,  consider 
a  piece  of  some  material  in  form  of  a  rod,  the  two  ends 
of  which  are  kept  at  different  temperatures  and  the  sides 
are  insulated.  Neglecting  dependance  of  VJU  on  tempera¬ 
ture  and  assuming  all  forces  except  VT  and  V<rf  are  zero, 
we  get: 

“  FO-’Jl^f+Wy  2„)  g  ^2 

-  Wl  7^  "  F  (U43^f  *  Iff  2 . 13 


and  current 

1  =  +  ***i)^f'  ”  +  ]  2.14 
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when  there  is  no  current  flow,  we  get 


L4|  L41 ; 


L-  js  i p  +  l— q  »y^2h  +  H.^2, 


2.15 


this  is  the  thermoelectric  power  describing  the  potential 
per  unit  of  temperature  difference  under  the  condition  of 
zero  electrical  current. 

Many  o.f  the  coefficients  L,\"  have  been  measured  in 
other  studies.  For  instance,  the  solvent-electric  current 
flow  coupling  coefficients  are  well  known  from  studies  by 
geophysicists  of  streaming  potential  and  electro-osmosis 
because  of  their  effects  on  self -potentials  in  drill-holes 
(Wyllie,  1955). 

The  coupling  effects  between  solute  flow  and  electric 
current'  known  as  diffusion  potential  are  also  well  known 
because  of  their  great  importance  in  electrochemistry. 

The  thermoelectric  properties  are  less  well  known,  but  a 
great  effort  is  being  made  to  study  them  because  of  the 
possible  application  in  direct  energy  conversion,.  Many 
measurements  of  electrokinetic  and  thermoelectric  coeffi¬ 
cients  have  been  made  by  Madden  and  his  group  (Kermabon, 
1956;  Madden  and  Marshall,  1959;  Madden  and  Cantwell,  1962). 

Section  III  -  Discussion  of  the  Coefficients  L,j 

In  this  section  we  discuss  the  physical  significance 
of  some  of  the  phenomenological  coefficients,  and  their 
relationship  to  the  more  familiar  electrochemical  and 
thermal  quantities.. 

Section  Ilia  -  Mobility 


The  mobility  U°  is  defined  as  the  speed  of  an  ion 


20 


in  infinite  dilution  under  a  unit  electric  field  when  all 
other  forces  are  zero.  The  velocity  of  a  positive  ion 
is  given  by3y/£p  and  that  of  a  negative  ion  is  given  by-Sn/c* 
where  and  3*  are  in  moles/cm^-sec,  and  Cp  and  C„ 
are  concentrations  in  moles/cm^.  The  minus  sign  for  nega¬ 
tive  ions  is  necessary  since  they  move  in  the  oposite 
direction  under  the  same  electric  field.  Prom  equations 
2.11,  we  find 


'cp  V=-(  Cp  2,16 

U=  (-S A  -F(L  id 

£  2.17 

In  practice  for  dilute  solutions,  the  coupling  coefficients 
and  are  small  compared  to  coefficients  \->>5  and 
,  and  we  have  therefore 


U  2  FL- aZf 


"IT  tfc  —  F 
*  r* 

>—  Vl 
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In  Table  2.1  typical  values  of  ionic  mobility  at  25°C 
at  infinite  dilution  are  listed  (Glasstone,  1954;  p.60): 


H 


r+ 


K 


.+ 


B 


,+■ 


N 


t+ 


Table  2.1  -  Ionic  Mobility  u  at  25°C 


36.2  X  10"H  OH" 

7.6  X  10-4  .  S0^“ 

6.6  X  10"4  CP" 

5.2  X  10-4  N0^ 

4.0  X  10"4  HCO^ 


(cm/sec) 

20.5  X  10-4 

8.3  x  10"4 

7.9  X  10"4 

7.4  X  10"4 
4.6  X  10"4 
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These  values  of  ionic  mobilities  indicate  that, 
apart  from  hydrogen  and  hydroxyl  ions,  most  other  ions 
have  comparable  velocities  around  4  to  8  X  10"^  cm/sec., 
therefore,  one  would  expect  that  hydrogen  and  hydroxyl 
ions  may  contribute  heavily  to  the  current  flow  and 
electrochemical  coupling.  But  as  we  shall  see  in  Section 
IVd  of  this  chapter,  when  several  ionic  flows  are  taking 
place  simultaneously,  the  contribution  of  each  ion  to 
the  current  flow  depends  on  its  concentration  gradient 
as  well  as  its  mobility.  In  general,  the  concentrations 
of  hydrogen  and  hydroxyl  ions  are  very  low  compared  to 
other  ions  such  asC*^  K+,  #l  Ma^and  therefore,  contribute 
little  to  actual  current  flow. 


Section  Illb  -  Electrical  Conductivity 

Electrical  conductivity  or  is  the  conductivity 
due  to  an  electric  field  ~V(p  when  all  other  forces  are 
zero,  and  it  is  given  by  equation  2.10..  Inserting 
and  3*  from  2.11  into  2.10  and  using  Ohm's  law 


T  =- 


<5T  V 


2.19 


we  get 


-  f(WUp'^C„0h)  2.20 


Section  Illo  -  Thermal  Conductivity 

In  the  absence  of  all  forces  except  for  temperature 
gradient,  the  flow  of  heat  is  governed  by  Fourier's  law 

p 

which  states  that  heat  flow  in  cal/cm  -  •  is  proportional 


if 
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to  the  temperature  gradient 

3fo  =  -  K  VT  2.21 

Comparing  equation  2.21  with  equation  2.11,  we  see  that 
for  v  P=  'u  =  t>  we  have 


but  in  general  the  chemical  potential  jUl  depends  strongly 
on  temperature  and  therefore  the  relationship  is  more 
complicated.  Taking  into  account  the  dependance  of  U  on 
temperature  and  ignoring  for  the  present  time  its 
dependance  on  pressure  and  concentration,  we  can  write 

VyU;  £  ^  VT  -  -  TS’  2Z  2.23 

where  is  the  partial  molal  entropy.  Inserting  2.23 
in  2.11,  we  get  for  heat  flow 

•=  —  —  2.24 

Therefore 

K  =  +  (^  -  L135p  -LwSh)  2.25 

Even  though  the  thermal  force  in  Onsager's  equation  is 
-vt/t  ,  we  have  assumed  tacitly  that  the  coefficient  K 
rather  than  L)(  is  a  constant.  This  is  an  agreement 
with  the  empirical  equation  2.1.  Actually,  this  assump¬ 
tion  is  not  always  true  and  for  some  material  the  coeffi¬ 
cient  K  decreases  with  increasing  temperature  which 


would  indicate  L(|  is  more  nearly  constant,  but  anyhow 
in  our  arguments,  none  of  these  coefficients  are  required 
to  be  strictly  constant. 

In  Table  2.2,  the  thermal  properties  of  some  common 
substances  are  listed  (Carslaw  and  Jaeger,  1947;  p.382) . 
Thermal  conductivity  K  is  given  in  calories/cm-sec-°C 
and  thermal  diffusivity  Si  (=  where  y?  is  the 

density  and  C  is  the  specific  heat)  is  given  in  cm  /sec. 

Table  2.2  -  Thermal  Properties  of  Some  Common  Substance 

(C.G.S.  Units.) 


Substance 

Conductivity  K 

Diffusivity 

Copper 

0.93 

1.14 

Air 

0.000058 

O.I87 

Granite 

0.006 

0.011 

Limestone 

0.004 

0.007 

Sandstone 

0.006 

0.011 

Soil  (average) 

0.0023 

0.0046 

Water 

0.00144 

0.00144 

Thermal  diffusivity  is  the  parameter  on  which  the 

rise  in  temperature'  depends..  This  is  so  because  a 

material  must  not  only  have  a  high  thermal  conductivity, 

but  it  must  also  have  a  very  low  specific  heat  coefficient. 

Otherwise,  much  of  the  heat  that  is  traversing  a  unit 

volume  will  be  used  to  raise  its  temperature  and  very 

little  will  be  passed  to  the  next  unit  of  volume.  The 

values  shown  in  Table  2.2  indicate  that  the  diffusivity  in 

rocks  is  very  small  and  therefore  the  heat  distribution  in 

rocks  is  quasi-static  when  our  length  dimensions  are  of 
4  q 

the  order  of  .10  or  Kr  cm.  A  given  picture  of  heat 
distributions  is  expected  to  remain  essentially  the  same 
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over  a  period  of  days,  months  or  even  years  as  we  already 
have  shown  in  Chapter  I.  For  instance,  in  the  case  of  an 
underground  nuclear  blast,  the  zone  of  high  temperature 
which  might  be  150  feet  in  radius  will  not  change  by 
more  than  a  few  per  cent  over  a  period  of  a  year  or  two. 
This  point  is  essential  to  the  usefulness  of  the  arguments 
presented  in  this  thesis,  otherwise,  the  effects  would 
be  so  short  lasting  that  they  would  be  of  no  practical 
application . 

For  other  equations  relating  diffusion  coefficient, 
transport  number,  and  Einstein's  relation  to  phenomeno¬ 
logical  coefficients  Lqv  ,  the  reader  is  referred  to 
Madden  (Madden  and  al.,  1957>  p.24-27). 

Section  IV  -  Application  of  Irreversible  Thermodynamics 
i  to  Study  Coupling  Effects  in  the  Earth 

Equation  2.11  is  used  in  this  section  to  derive 
specific  equations  for  coupling  of  a  pair  of  flows. 
Strictly  speaking,  it  is  impossible  to  isolate  a  pair 
of  flows  and  their  corresponding  forces,  because  the 
presence  of  a  gradient  X  will  cause,  through  coupling, 
the  establishment  of  other  potential  gradients  and  flows 
as  well.  Nevertheless,  this  assumption  is  believed  to 
be  a  good  one,  because  it  allows  us  to  isolate  one  effect 
at  a  time  and  study  it  in  detail.  Then  for  the  case  of 
several  couplings,  it  is  easy  to  superimpose  the  effects 
since  equation  2.11  is  linear. 

In  one-dimensional  flow  studies,  the  formal  solu¬ 
tion-  of  the  coupling  equations  is  straight  forward,  and 
it  is  further  simplified  because  one  of  the  flows  is 
usually  assumed  to  be  zero.  In  studying  coupling  effects 
in  the  earth,  we  are,  in  general,  forced  to  work  in  three 
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dimensions,,  and  furthermore,  we  are  often  interested  in 
transient  conditions.  In  what  follows,  we  shall  assume 
that  the  time  scale  for  electric  current  is  much  shorter 
than  the.  time  scale  for  diffusion,  heat,  and  solvent 
flow,  so  that  the  charge  distribution  is  set  up  instan¬ 
taneously.,  With  this  assumption,  we  have  in  the  absence 
of  electrical  sources 

7.1=0  2 .26 


Section  iva  -  Electrokinetlc  Coupling 

To  investigate  the  effect  of  electrokinetlc  coupling 
on  currents  generated  in  the  earth,  we  shall  assume 
all  forces  except  the  driving  force  -VP  and  the  induced 
force  -V  are  zero.  The  dependance  of  on  pressure 
can  be  neglected  (Marshall,  1959;  P*33).  Combining  Jp 
arid  to  give  the  current  flow  J  ,  we  get 


=  -WlvP  -  L*,  V0  2.27 

I  -  -  L***  VP  -  Uv  2  28 


where 


W3  -  2r  •+• 

2.29 


Here  we  have  used  Onsager's  relations  2.9.  The  meaning 
of  the  coefficients  in  equations  2.27  and  2.28  are  given 
here. 


electrical  conductivity  0“ 
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=  permeability  K 

w  p 

=  streaming  potential  coefficient 

O 

=  electro-osmotic  coefficient 
o 

Onsager's  relation  allows  us  to  get  both  stream¬ 

ing  potential  and  electro-osmotic  coefficients  from  one 
single  .measurement  of  coupling  effect. 

Experimental  data  on  streaming  potential  coeffi¬ 
cients  are  scarce.,  but  recently  Madden  and  his  group  have 
made  systematic  measurements  of  streaming  potential  on 
rock  samples  (Kermabon,  1957;  Madden  and  Marshall.,  1959; 
Madden  and  Cantwell,  1962).  In  Table  2.3  are  listed  the 
summary  of  the  measurements  made  on  rocks  by  Kermabon, 
and  Madden,  and  on  clays  by  Olsen.  These  measurements 
were  made  using  tap  water.  At  the  present  time  these 
studies  are  being  continued  to  determine  the  effects  of 
variations  in  salinity  on  the  streaming  potential  in 
geologic  materials. 

From  the  magnitudes  of  the  streaming  potential 
coefficients  in  Table  2.3,  we  see  that  electrokinetic 
couplings  can  be  very  important  in  the  earth.  Streaming 
potential  coefficients  for  geologic  materials  are  as 
high  as  10  dr  20  mv/atm.,  and  under  favorable  conditions 
even  a  pressure  difference  of  only  a  few  atmospheres  can 
produce  detectable  differences  of  potentials  of  a  few 
hundred  millivolts.  These  favorable  conditions  depend 
largely  on  the  geometry  which  will  be  discussed  in  the 
following  chapter. 
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The  continuity  equation  for  incompressible  solvent 
flow  is 

V.  Js  =  5 (r,t)  2.30 

This  equation,  which  is  valid  for  uncoupled  flows,  remains 
valid  when  couplings  are  present  as  long  as  3$  ,  repre¬ 

sents  the  total  contributions  due  to  all  thermodynamic 
forces  (DeGroot,  1961;  p.95). 

Equations  2.26,  2.27,  2.28,  2.30  and  the  boundary 
conditions  on  the  continuity  of  potentials  and  normal 
components  of  flows  are  sufficient  to  solve  for  any  geom¬ 
etry  by  standard  techniques  used  in  solving  boundary-value 
problems.  But  this  would  cause  unnecessary  difficulties 
because  we  have  two  coupled  partial  differential  equations 
with  coupled  boundary  conditions.  A  great  simplification 
can  be  obtained  if  we  remember  that  vP  is  the  cause  and 
is  the  effect,  and  while  we  have  to  keep  both  terms 
on  the  right  hand  side  of  the  equation  2.28,  we  can 
neglect  the  second  term  on  the  right  hand  side  of  equation 
2.27  ( V  P  »  Vj/  )•  A  quick  calculation  will  convince 
us  of  the  validity  of  this  assumption.  Taking  some 
typical  values  from  Table  2.3,  we  have 


L22 

=  10TX° 

(MKS) 

l 

l23 

\ 

=  L32  1 

=  10~10  (MKS) 

i 

L33 

-  :..io"3 

(MKS) 

AT 

=  5  atm 

=  5  X  105  (MKS) 

=  50  mv 

»  5  X  10"2volts  (MKS) 

therefore 
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Lj' =  S  xid~'d 

l,;  hfi  =  sx\o,z 

L &P  =>  IC>* 
u,;  Ajz (  =  ** ID_ 


Prom  these  calculations,  it  is  clear  that  A^  is 
much  smaller  than  AP  while  Lj^AP  is  comparable 
to  L  «  A  ^  .  This  amounts  to  uncoupling  the  solvent 

flow  equation.  We  then  get 

3*  Qt  -U^VP  2.31 

1  »  -  L»%*  Vf  2 .32 

where  -f  is  defined  as  the  sum  of  electric  potential 
and  the  potential  induced  through  coupling 

f  =  $  +  (—A.)  P  2.33 

'  Li3J  ' 

and  using  2.26,  we  get  the  equation  satisfied  by  £ 

Ve  f  =  O  2.34 

Equation  2.31  for  the  solvent  flow  is  separated  and 
can  be  solved  by  straight  forward  methods.  Equation  2.32 
gives  the  current  flow  in  terms  of  a  potential  f  which 
we  shall  call  the  "total  electric  potential"  and  is  made 
up  of  the  ordinary  potential  (f)  and  a  contribution  due 
to  coupling  (L^/L,; )  f  •  Potential  f  satisfies  the 
Laplace's  equation  in  each  medium,  but  its  value  on  the 
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boundary  depends  on  y  and  P  ,  which  amounts  to  saying 
that  £  Is  caused  by  sources  distributed  along  the  bounda¬ 
ries  Equation  2.33  shows  clearly  the  dependance  of  f  on 
the  streaming  potential  coefficient  (  )  . 

In  the  next  two  sections  of  this  chapter,  similar 
equations  are  derived  for  the  thermoelectric  and  electro¬ 
chemical  coupling  cases  and  are  shown  to  be  identical  to 
equations  2.32,  2.33,  and  2.34. 

Section  IVb  -  Thermo-Electric  Coupling 

A  similar  analysis  is  used  to  investigate  the  thermo¬ 
electric  coupling  effect,  but  as  we  pointed  out  when 
discussing  thermal  conductivity  K  ,  the  thermoelectric 
coupling  is  slightly  more  complicated  because  of  the  role 
played  by  the  chemical  potentials.  Using  equations  2.11 
and  2.23,  we  can  write 

3^  =  -  L.'  VT  2.35 

X  =  -L  2.36 

_  L,y  Sv,j 

2  37 

F  [b L+h>A.  _  (Lis?r<L„?.)Sf  -  (L 

C"  = 

,  / 

In  this  case  the  coefficients  L()  and  U*,  are  not  equal, 
and,  therefore,  one  single  measurement  of  coupling  effect 
is  not  enough  to  give  us  both  coefficients,  but  rather 


where 


u:  = 


14  * 

/ 

53 

L3>  a 
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two  distinct  measurements  of  thermoelectric  e.m.f .  and 
Peltier  heat  are  necessary 


l;  , 

fc&L  ■ 

thermal  conductivity  K 

/ 

/Vj»=o 

n 

0 

(i 

_ H 

''aiK 

thermoelectric  e.m.f.  9 

Peltier  heat  Tf 

^ 1  4t=. 

In  Table  2.4  are  listed  some  of  the  measured  thermo¬ 
electric  coefficients  for  geologic  materials  (Madden,  1962) . 
These  measurements  were  made  by  comparing  the  voltage 
developed  across  the  rock  sample  when  a  temperature  gradient 
exists,  to  that  developed  across  a  tube  of  water  to  elimi¬ 
nate  the  temperature  effect  on  the  recording  electrodes. 
Therefore,  these  values  represent  only  the  difference  in 
the  coupling  coefficients  of  the  rock  and  water,  and  the 
actual  coefficients  cannot  be  known  without  a  knowledge  of 
the  thermoelectric  properties  of  water.  But  as  we  shall 
see  later  on,  we  are  only  interested  in  the  differences  in 
coupling  properties . 

Table  2.4  -  Measured  Thermoelectric  Coupling  Differences 


Rock  Type 

Code  Sample* 

Thermoelectric 

coefficient** 

Y4 

.257 

Y7B 

.475  ■ 

.  Y9 

.002 

Sedimentary  Rocks 

Yl8 

.152 

Colorado  Plateau 

Y3^ 

.107 

Y36 

.135 

Y4l 

.068  ' 

Y49 

.228 

Y52 

.158 
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SA  7 

.195 

SA  8 

.196 

SA  9 

.348 

Altered  Volcanics 

SA  14 

.120 

Arizona 

SA  27 

1.36 

SA  29 

.070 

SA  31 

.076 

SA  37 

.350 

SA  40 

.449 

DC  3 

.268 

BC  4 

.150 

Quartz -Latite 

BC  5 

.444 

Porphyry,  New  Mexico 

BC  9 

.375 

BC  10 

.184 

S  2 

-.08 6 

S  3 

.108 

Dakota  Sandstone 

S  7 

.081 

S  10 

.243 

S  21 

1.12 

♦Samples  codes  are  those  used  by  T.  R.  Madden 

♦♦Values  listed  represent  the  difference  between  the  thermo- 

electrlc  coefficients  of  rock  samples  and  water, _ 

From  the  data  in  the  above  table,  it  is  obvious  that 

the  thermoelectric  effect  can  be  very  important.  A 

temperature  difference  of  100°C  can  produce  sizeable 

anomalies  of  20  to  100  millivolts.  This  explains  why  at 

the  present  time  so  much  effort  is  spent  in  studying  the 

thermoelectric  effect  for  direct  energy  conversion. 

However,  for  the  case  of  an  underground  nuclear  blast, 

this  effect  is  believed  to  be  less  important  than  the 

streaming  potential  effect  because  of  the  even  larger 

electrokinetlc  coupling  coefficients. 

For  similar  reasons  to  those  given  in  Section  IVa, 

the  continuity  equation  for  heat  flow 


v-3®  *  /c1f  - 


2.38 
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remains  valid  when  couplings  are  present.  Furthermore, 
making  the  assumption  Ltl' V T  when  temperature  T 

is  the  driving  force,  we  get 

Ja  a  _  l'  vt 
X  =  -  -vf 

where  T  Is  related  to  T  and  by  the  relation 

*  -  ^  +  (i,)T  2A1 

and  satisfies  the  Laplace's  equation 

^*“4  =  0  2-42 

Clearly,  equations  2.40  and  2.42,  giving  the  current 
flow  when  thermoelectric  couplings  are  present,  are 
identical  with  equations  2.32  and  2.34  which  give  the 
current  flow  when  electrokinetic  couplings  are  present. 

A  quick  calculation  will  convince  us  that  the  assumption 
L,'  VT »  L-ii  is  a  good  one.  Taking  some  typical  values 

L>-  ^  ^  /^\-yOe  c-  °C 

=  W.T  =  I.7S-  °c/i*U 

&T  =  80  °C  =  O.  I  \jo^- 

we  get 

l.[h>T=Zoo  O.V.J  =  .ns- 

This  last  assumption  is  very  helpful  because  it  decouples 
the  equations  2.38  and  2.39  for  heat  flow  (but  not  the 


2.39 

2.40 
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equations  2.40  and  2.42  for  current  flow,  because  f 
depends  on  both  and  T  .),  and  therefore,  we  can 
solve  the  heat  flow  problem  by  standard  methods  and  use  the 
results  to  obtain  the  induced  current  flow. 

Section  IVc  -  Electrochemical  Coupling 

In  both  electrokinetic  and  thermoelectric  couplings, 
the  flow  of  current  depended  on  cross-coupling  terms 
and  L/  ,  but  it  is  a  quite  different  matter  when  diffusion 
of  ions  is  taking  place.  This  is  because  the  ions  them¬ 
selves  carry  the  electric  current,  so  that  the  coupling 
is  direct,  concentration  gradients  and  electrical  poten¬ 
tial  gradients  appear  together  as  the.  driving  force.  We 
then  have 


~  ^*0  v/ 

2.43 

-  —  *~ifS  ~ 

-  LwZ,,)  Vfi 

2.44 

As  stated  above  the  direct  coupling  effects  are  so 
important  that  the  cross-terms  can  be  neglected.  We 
further  assume  ,  since  the  concentrations  of  the 

positive  and  the  negative  ions  can  not  be  very  different, 
otherwise,  tremendous  electric  fields  would  be  set  up 
(this  ils  true  only  when  two  kinds  of  ion,  one  positive  and 
one  negative  are  present) .  (The  more  general  case  is 
discussed  later  on.)  With  these  assumptions,  we  have 


I 


2.45 
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where  the  chemical  potential  JLL  is 

-Ac  =  J^oL  +  CX  2.46 

Using  equation  2.26,  we  get 

i  =  _  L,;vf  2,47 

where 

V  =  </  -t  a  £.48 

and  t 

V2  f  -  o  2.49 

The  above  equations  are  in  essence  identical  to 
equations  governing  electrokinetic  and  thermoelectric 
couplings.  The  coefficient  multiplying  ju.  in  equation 
2.48  can  be  considered  as  the  electrochemical  diffusion 
coupling  coefficient. 

It  is  remarkable  that  the  equations  governing  all 
three  effects  are  identical  and  therefore  all  three 
effects  can  be  studied  together.  The  actual  difference 
between  the  three  effects  is  in  the  physics  of  the  problem 
which  results  in  different  kinds  of  boundary  conditions. 

Section  IVd  -  Several  Ionic  Flows 

We  can  still  neglect  the  cross-coupling  terms  when 
several  ionic  species  are  present.  Using  2.46,  we  have 

FL;^v/  2.50 

2.51 

2.52 


1+  =  - 
1/JL;*  KT  JC,- 
L;  2  Ui  Cf-  * 

Fi li 


-  37  - 


where 


0 

*-  =  ion  of  the  ith  kind  -?/=  charge  of  ith  ion 

Ut*  =  mobility  C**»  concentration 

^*i’  —  (i  =  3*4,...,n) 

Current  flow  T  is  given  by 


1  -  -  F  Z  Wii  -  F8(IL;^J  v/ 


2.53 


Using  2.51  and  2.52,  we  get 


1=  -  F^IW)  V/4T 


I  u;  y^r  ' 
ru.c.-z,-  . 


2.54 


which  indicates  that  the  contribution  to  the  current  flow 
by  each)  ion  is  a  "weighted"  average  of  its  mobility  and 
concentration  gradient.  The  ion  with  the  greatest  mobility 
and  concentration  gradient  contributes  most  of  the  current 
flow.  •; 


Section;  IVe  -  Summary  of  Coupling  Equations 
i  ' 

Starting  from  Onsager's  relationships,  we  were  able 
to  derive  a  generalized  equation  for  the  current  flow  due 
to  electrokinetic,  thermoelectric,  or  electrochemical 
coupling.  For  each  medium  we  have: 


2.55 


V‘f  =  o 

and  the  electric  potential  p  Is  given  by 

^  i -  CY 


2.56 


2.57 
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Using  the  boundary  conditions  for  the  continuity  of 
electric  potential  ( i-A  )  and  normal  component  of 
current  flow  (3,.n  we  can  write  the  boundary 

conditions  for  f : 


2.58 


and 


Cr  zjj 


*  "isn 


2.59 


where  subscripts  1  and  2  refer  to  the  two  media  in  contact; 

~Y  i;s  the  generalized  driving  potential  (  ,  T  or ); 

Xi  is  the  value  of  this;  potential  at  the  boundary 
(  P  y  T  and  JUL  are  continuous  at  the  boundary);  -f  is 
the  "total  electric  potential"  giving  rise  to  current  flows 
(  **  )  is  the  conductivity;  and  C  is  the  ratio  of 

the  cross-term  to  electrical  conductivity.  We  shall  call 
C  the  "coupling  coefficient." 

Potential  f  is  uniquely  determined  by  the  equation 
2.56  and  the  two  boundary  conditions  2.58  and  2.59.  It  is 
important  to  note  that  the  driving  potential  y  enters 
only  in  the  boundary  condition  2.58  through  its  value  at 
the  boundary  Yyl  ,  and  its  value  inside  each  medium  is 
Immaterial.  In  other  words,  the  generation  of  currents 
depends  essentially  on  interaction  between  the  equipoten- 
tial  surfaces  and  the  boundaries.  This  interaction  creates 
a  series  of  "sources"  at  the  boundaries.  Therefore,  we  can 
expect  ^that  the  geometry  of  the  coupling  differences  will 
play  a  dominant  part  in  the  magnitude  and  the  distribution 
of  the  observed  effects,  while  other  factors  such  as  the 
geometrical  distribution  of  conductivities  will  influence 
much  less  the  electric  potential  picture. 


CHAPTER  III 


THEORETICAL  SOLUTIONS  FOR  SOME  SIMPLE  GEOMETRIES 

In  the  previous  chapter,  using  the  principles  of 
the  thermodynamics  of  irreversible  processes,  we  derived 
the  equations  for  the  potential  distribution  due  to  the 
coupling  of  flows.  In  this  chapter,  these  equations  are 
solved  analytically  in  terms  of  the  driving  potential 
(or  rather,  its  value  at  the  boundaries)  for  the  particu¬ 
larly  simple  geometries  of  horizontal  and  vertical 
layering.  The  problem  of  random  distribution  of  inhomo¬ 
geneities  is  solved  for  one -dimensional  variations . 
Numerical  computations  are  carried  out  for  concentrated 
heat  ahd  pressure  sources  buried  in  a  horizontally 
layered  half  space;  concentrated  heat  sources  in  the 
presence  of  a  buried  slab;  and  depth  dependent  chemical 
potentials  in  the  presence  of  vertical  contacts  and 
dikes. . 

The  laboratory  measurements  discussed  in  the  previous 
chapter  show  how  efficiently  rocks  can  couple  pressure, 
temperature,  and  chemical  potential  gradients  into  elec¬ 
tric  potential.  At  the  earth's  surface,  it  can  be 
assumed  that  these  gradients  do  not  exist  or  are  neglible, 
either  because  the  air  acts  as  a  short  circuit  for  these 
gradients  (for  example,  air  has  an  essentially  infinite 
permeability,  and  as  a  result,  the  fluid  pressure  at  the 
surface  of  the  earth  will  be  the  same  everywhere),  or  else 
because  the  flowB  are  confined  to  the  near  vicinity  of 
the  underground  source.  However,  surface  currents  are 
possible  if  there  is  coupling  at  depth  as  will  be  shown 
shortly. 
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Section  I  -  Horizontally  Layered  Earth 

To  begin  with,  we  consider  the  Idealized  but  realis¬ 
tic  case  of  a  layered  earth.  Our  model  (Pig.  3*1)  consists 
of  homogeneous,  conductive  layers  overlying  a  homogeneous 
half  space.  Each  layer  has  a  thickness  ,  electrical 
conductivity  O".  ,  and  coupling  coefficient  C-  .  For 

such  a  layered  earth  composed  of  n,  layers,  we  have: 


=  °  (  *  *  »/*,•••,«) 


3.1 


and  boundary  conditions 


\>i 


3.2 


V,U»-| 


3.3 


Because  of  the  axial  symmetry,  it  is  convenient  in  this 
case-  tp  use  cylindrical  coordinates  )  with  the 

£  -axis  pointing  downward.  We  shall  treat  in  detail 
three  cases — homogeneous  half  space,  one  layer  over  a 
half  space,  and  two  layers  over  a  half  space — to  give  a 
clear  idea  of  the  anomalies  to  be  expected  at  the  surface. 

Homogeneous  Half  Space  -  Air  is  a  very  poor  conductor  of 
electricity  compared  to  the  earth,  and  therefore  the 
conduction  currents  in  the  air  are  negligible.  The  only 
boundary  condition  that  the  potential  f  must  satisfy 
is  therefore: 


*1  2-0 


in 


o 


3.4 


0 


FIGURE  3.1  HORIZONTALLY  LAYERED 
HALF  SPACE 


0 


FIGURE  3.2  ONE  LAYER  OVER  HALF  SPACE 
GEOMETRY  FOR  THE  THERMO¬ 
ELECTRIC  CASE 
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-f  satisfies  Laplace's  equation,  and  by  3.4  its 
normal  derivative  must  vanish  on  the  boundary.  One 
possible  solution  satisfying  these  two  conditions  is 
-f  =  constant,  and  using  the  uniqueness  of  solutions 
theorem  (Ramsey,  1959;  p.H7>.  problem  2) ,  we  see  that 
f  must  be. constant  everywhere  in  the  half-space; 
therefore,  there  are  no  currents  generated.  The  elec¬ 
tric  potential  which  is  what  we  try  to  measure  (see 
Appendix  B)  is  given  by 

<f>  =•  f  -  c  y  3.5 

where  Y  is  essentially  constant  on  the  surface.  For 
instance  if  ~Y  represents  .temperature  rise,  it  is  zero  on 
the  surface  because  the  rocks  are  poor  conductors  of 
heat,  and  as  a  consequence,  the  zone  of  high  temperature 
is  confined  to  the  near  vicinity  of  the  underground 
source.;  Even  if  the  temperature  rise  reached  the  surface 
air  which  is  a  very  good  conductor  of  heat  (compared  to 
the  roqks)  would  act  as  a  short  circuit  to  bring  the 
temperature  everywhere  on  the  surface  to  its  background 
value.  If  the  driving  force  is  the  pressure  1?  ,  it 

will  be  zero  everywhere  on  the  surface  of  the  ground 
because  air  has  infinite  permeability.  The  chemical 
potential is  also  constant  on  the  surface,  or  we 
hope  so,  because  its  value  depends  on  temperature,  pres¬ 
sure  and  concentration  of  ions  which  depends  on  the 
oxidation  conditions  at  the  surface.  As  long  as  one  is 
making  measurements  in  the  same  medium  the  chemical  poten 
tial^^  should  not  change  laterally.  Therefore,  accord¬ 
ing  to  3.5,  the  electric  potential  <j/>  is  the  same 
everywhere,  and  no  variations  can  be  detected. 
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The  above  negative  result  is  of  tremendous  importance, 
because  it  Indicates  that  even  though  there  is  a  driving 
force  in  the  medium,  and  there  is  coupling  between  the 
flows  as  predicted  by  Onsager's  relations,  the  homogeneity 
of  the  medium  cancels  on  the  surface  any  effect  that  may 
be  produced  at  depth.  This  is  because  there  is  no 
boundary  cutting  across  the  equipotential  surfaces  and 
inducing  electrical  sources.  The  following  cases  will 
show  that  effects  detectable  at  the  surface  arise  whenever 
inhomogenelties  in  rock  properties  cut  across  the  equi¬ 
potential  surfaces  associated  with  the  various  driving 
forces . 

Incidentally,  equation  3*5  implies  that  as  long  as 
both  potential  electrodes  used  in  measuring  $6  are  in  the 
same  mddium,  it  is  impossible  to  distinguish  between  the 
potential  difference  and  the  "total  potential"  differ¬ 

ence  A-f  (we  are  assuming  A/  —  o  on  the  surface).  The 
two  potential  differences  are  distinct  whenever  the  elec¬ 
trodes  are  in  different  media.  We  then  have 

ap/  =  Af-  my  3.6 

where  AC  is  the  difference  in  coupling  coefficients. 

One  Layer  Over  a  Half  Space  -  For  the  case  of  a  layer  of 
thickness  d  overlying  a  half  space,  we  have  the  follow¬ 
ing  boundary  conditions 
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Using  cylindrical  coordinates  (Morse  and  Feshbach, 
1953) ,  the  solution  can  be  written  in  the  following  form 


f  =  3.9 

o 

and  using  the  boundary  conditions  3.8  and  3.7,  we  get  (see 
Appendix  A) 


ft 


3.10 


where  is  the  electrical  reflection  coefficient,  so 
familiar  to  people  using  resistivity  methods,  and  )7(A) 
is  the'  Fourier-Bessel  transform  of  the  driving  potential 
"yCr,  i)  evaluated  at  o(. 


K 


*-\ 


<5y-ST 

(T,  *«-. 


3.11 


Xx  *  j  YM)  * *r  3 . 12 

o 

Therefore,  we  get  for  the  potential  -f-  in  Medium  I 


3.13 


where 

A-  3.14 

Clearly  enough  the  value  of  £  (vj  JZ)  in  equation  3.13 
depends  only  on  the  value  of  the  driving  potential 
evaluated  at  the  boundary  between  the  layer  and  the 
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conductive  half  space.  The  effect  of  denominator  A.  is 
to  give  rise  to  multiple  reflections  (series  of  images) 
due  to  contrast  in  conductivity  between  the  two  media.  If 
the  resistivity  contrast  is  zero  (  -  o  ,  L*  I  ),  the 

equation  3.13  giving  f  is  slightly  simplified  but 
remains  essentially  unchanged. 

Assume  that  the  boundary  between  the  layer  and  the 
half  space  does  not  cut  across  any  equipotential  surface 
connected  to  the  driving  force.  In  this  case 

=  x 

where  Y  is  a  constant.  Using  equation  3.12,  we 
have 

-  X  [  ^(Xr)T<Jr  -=  3.15 

Substituting  from  3.15  into  3.13>  we  get: 

■fr,(r,2)  =  C c.-O  X  3.16 

and  A-f,  =  O 

This  is  in  agreement  with  what  was  stated  while  discussing 
the  case  of  a  homogeneous  half  space.  If  the  boundary 
does  not  cut  across  equipotential  surfaces,  as  it  is  the 
case  here,  no  effect  can  be  observed  on  the  ground.  This 
result  has  nothing  to  do  with  our  choice  of  horizontal 
layering. 

Application  to  Buried  Heat  Source  -  One  of  the  long  last¬ 
ing  aftereffects  of  an  underground  nuclear  blast  is  the 
release  of  a  large  amount  of  heat,  and  therefore  creation 
of  a  temperature  gradient.  The  heat  released  is  roughly 
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5  X  10 ■l’2  joules/kiloton  of  TNT  (Johnson  and  al.,  1959), 
and  since  the  rocks  are  poor  conductors  of  heat,  this 
energy  will  be  confined  to  the  near  vicinity  of  the  shot 
for  an  extremely  long  period  of  time  if  aeration  and 
ground  water  movements  are  negligible.  Most  of  the 
energy  goes  into  raising  the  temperature  of  the  water 
present  to  its  boiling  point  and  turning  some  of  it  into 
steam.  This  steam  will  be  in  equilibrium  with  water, 
and  therefore,  its  temperature  will  be  the  boiling  point 
temperature  of  water  at  the  prevailing  pressure.  In  this 
manner  no  part  of  the  high  temperature  zone  will  be 
above  the  boiling  point  of  water  (Johnson  and  al.,  1959) j 
but  rather  there  will  be  a  roughly  spherical  zone  at  a 
uniformly  high  temperature  and  a  surrounding  medium  at 
ambient  temperature.  The  radius  of  the  high  temperature 
zone  for  20  KT  blast  can  be  as  much  as  150  to  200  ft. 

In  our  model  (Pig.  3*2)  the  spherical  heat  source 
at  temperature  is  centered  at  S  ,  and  a  is  the 

radius  , of  the  high  temperature  zone  cut  be  the  boundary 
at  .  The  background  temperature  is  assumed  to  be 

zero  degree,  because  assuming  a  background  temperature 
different  of  zero  will  not  change  anything,  as  only  the 
differences  in  temperature  are  important.  The  assumed 
temperature  distribution  along  the  boundary £sd  is  shown 
in  Fig.  3.3 
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With  this  idealized  picture  of  the  temperature  distri¬ 
bution  we  get: 


and 


%*,<>)  *  (***)&r  c0  T«  (  d)  5> 

*  n 


00  _  i 

-e 


0 


r.Cti)3.(f4)4  3.17 


In  going  from  3.13  to  3*17*  we  have  made  a  slight  change 
of  variable  for  computational  purposes. 

The  product  (c,-cz)Tp  in  the  above  equation 
represents  the  strength  of  the  induced  electrical  sources,  and 
depends  on  the  temperature  difference  T„  and  the  coupling 
difference  (  C,  -  C2  ) .  According  to  this  equation,  it  is 
possible  to  have  media  with  strong  coupling  properties  and 
still  observe  no  surface  effect,  if  the  coupling  coefficients 
are  equal.  An  asymmetry  of  the  coupling  properties  is 
essential  to  the  generation  of  electric  currents. 

The  ratio  ^ (r,o) /(<£,- (*)T0]  ,  which  we  shall  call  the 

"geometric  factor,"  is  a  dimensionless  number  and 
represents  the  effect  of  distance  X"  ,  thickness  of  the 
layer  cl  ,  size  of  the  shot  a  ,  and  the  conductivity 
contrast  Hz  1  •  A  high  speed  computer  has  been  used  to 
calculate  this  geometric  factor  for  several  values  of 

and  fi2| .  The  results  are  shown  in  Fig.  3»^« 

To  get  specific  answers  in  millivolts,  one  has  to  multiply 
the  ordinates. in  the  figure  by  the  appropriate  temperature 
rise  X  and  the  difference  of  coupling  (  c,-  Cz  ) .  It 
is  important  to  note  the  strong  dependence  of  the  geometric 
factor*  on  conductivity  contrast  as  well  as  the  thick¬ 

ness  d  of  the  overlying  layer.  This  is  because  the 
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thermoelectric  effect  is  highly  localized,  i.e.,  the  sources 
stop  abruptly  at  distance  cl  . 

Application  to  Buried  Pressure  Source  -  When  a  nuclear 
charge  is  detonated  underground  it  acts  as  a  pressure  source 
as  well  as  a  heat  source,  and  a  differential  pressure  of 
tens  or  even  hundreds  of  atmospheres  is  probably  the  case 
over  a  period  of  days  or  months .  We  now  apply  the  equa¬ 
tion  3.13  to  predict  the  instantaneous  potential  distribu¬ 
tion  due  to  a  concentrated  pressure  source. 

In  our  model  (Pig.  3*2)  the  spherical  pressure  source 
is  centered  at  S  at  a  depth  from  the  surface.  We 
shall  assume  that  the  earth  has  a  homogeneous  permeability. 
From  equations  2.30  and  2.31,  we  get: 

-  Li  V*F  -  s.V?-?.) 

and  remembering  that  the  fluid  pressure  must  be  zero  on  the 
surface,  we  get  the  solution 

P(r ,i  )  =  15-1  -  ILl 

R  ft' 

where  . 

R  -  l  i-vowi1'* 

and 

ft'  -  [rl+  (2+fc)m]4 

and  the  Bessel  transform  of  evaluated  at  the  bound¬ 

ary  f  sol  is 


Substituting  this  value  of  Cx)  in  equation  3*13,  we 
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get  on  the  surface  £-° 

r~  -H  -(.$**)$ 

o  \  -t-  5?.,  i* 


Expanding  the  denominator  and  integrating,  we  get 


V  rl4  (h+(ZH,t»)d)£-l 


3.16 


where  h  is  the  depth  to  the  pressure  center,  Jt  is  the 
thickness  of  the  layer,  r  is  the  distance  along  the 
surface,  is  the  conductivity  contrast,  and  a  is  the 

radius  at  which  the  pressure  is  measured  to  be  1?D  . 

Figure  3-5  shows  the  values  of  the  dimensionless  factor 

computed  for  several  values  of  (j),  (|)* 
and  K-i.i  .  Pressure  does  not  drop  off  as  abruptly  as 
temperature,  and  consequently  the  streaming  potential 
sourcesiare/qu^Lte  as  localized  as  the  thermoelectric  sources. 
This  is  indicated  by  the  slight  dependence  of  the  geometric 
factor  on  conductivity  contrast  /Ra,(see  Fig.  3*5). 

\ 

Two  Layers  Over  a  Half  Space  -  The  only  interest  in  this 
problem  is  to  show  what  happens  to  the  coupling  potentials 
when  more  than  one  boundary  cuts  across  the  equipotential 
surfaces.  See  Figure  3-6. 

Using  the  equations  and  the  boundary  conditions  3.1* 

3.2,  and  3.3*  we  get  (See  Appendix  A) 


too  .  (*«„ *6ryMI 


3.19 


where 


A.  =■ 


A  = 


FIGURE  3.4  PLOTS  OF  THERMOELECTRIC  POTENTIAL 


fi<r)/[(c,-c*)T0]  (SEE  FIG.  3.2  FOR 

EXPLANATIONS) 


and 
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Comparing  equation  3.19  and  3.13*  we  conclude  there¬ 
fore,  that  as  a  first  approximation  (l.e.,  )  the 

coupling  effect  in  this  case  can  be  considered  as  the  sum 
of  two  independent  effects  due  to  and  course, 

for  sources  of  equal  strength,  this  superposition  is 
weighted  by  the  distance  from  the  surface.  The  effect 
from  the  source  nearer  to  the  surface  is  the  more  important 
(this  is  suggested  by  the  terms  and  e  un^er 

the  Integral  sign) .  Consequently,  we  shallTestrict  our¬ 
selves  to  the  simple  cases,  because  we  can  get  a  rough 
idea  of  the  more  complicated  situations  by  superimpositing 
the  results  for  one  layer. 

Section  II  -  Buried  Slab 

.  -  ...  r  —  „  - — 

In  this  section  we  shall  compute  the  Induced  thermo¬ 
electric  potential  when  a  spherical  heat  source  interacts 
with  an  inhomogeneity  in  shape  of  a  buried  slab.  A  semi¬ 
infinite  slab  of  thickness  iw  Is  embedded  at  a  depth  -■  d. 
in  an. otherwise  homogeneous  half  space.  A  high  temperature 
spherical  zone  of  radius  R  is  located  at  the  same  depth 
at  a  distance  a  away  from  the  slab  (see  Fig.  3*7). 

Dipole  sources  are  induced  along  the  portion  of  the  bound¬ 
ary  which  is  intercepted  by  the  high  temperature  zone. 

Due  to  the  very  complicated  boundaries,  this  problem 
cannot  be  solved  analytically,  but  when  the  contrast  In 
conductivity  is  zero  (  =  o  ),  the  solution  can  be 

obtained  by  assuming  a  boundary  distribution  of  dipole 
sources  and  integrating  numerically.  Indeed,  for  , 

we  have: 
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3.39a 


Ofi  2b-— 


3.39b 


f,"  ^  v  "  Cc\ ~  XL 


3.39c 


According  to  these  equations  the  potential  -f  satisfies 
the  Laplace's  equation  in  each  medium  and  the  current  is 
continuous  everywhere,  while  the  potential  f  is  discon¬ 
tinuous  when  crossing  the  boundary.  This  means  that  the 
thermoelectric  potential  is  due  to  dipole  sources  distri¬ 
buted  along  the  boundary.  The  dipole  strength  is  (<,-£,)  7IZ 
per  unit  area.  The  dipole  potential  of  a  small  area  Afl 
located  at  (  %g  ),  and  having  the  orientation  h  is 


where 


1  ZTF 


3.40 


The  total  potential  -f,  is  obtained  by  summing  up  the 
contribution  to  all  the  dipoles.  This  problem  has  been 
solved  on  an  IBM  709  and  the  results  for  several  values 
of  vn 


ct ,  and 


d  are  shown  in  Fig.  3-10.  These  curves 


show  some  interesting  behaviours,  but  the  most  important 
feature  is  the  smallness  of  the  geometric  factor 
in  this  case  as  compared  with  the  case  treated  in  Section 
I  (Fig.  3.4). 


/ 


X 


FIGURE  3,7  A  BURIED  SLAB  AND  A  CONCENTRATED 
HEAT  SOURCE 


Oitfonct 

FIGURE  3.10  THERMOELECTRIC  POTENTIAL  fi  (»)  /  [(c  i-Ct)Toj0VER  A  BURIED  SLAB  FOR  R*,«0  (SEE  FIG  3.7 
FOR  EXPLANATION  I 
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Section  III  -  Vertically  Layered  Earth 

Another  idealized  but  useful  model  is  a  vertically 
layered  earth.  It  can  be  treated  analytically  and  it  is 
geologically  useful  as  a  model  for  vertical  contacts  and 
dikes  (see  Fig.  3*8  and  3.9). 

Geophysicists  have  observed  many  spontaneous 
polarization  anomalies  that  can  not  be  associated  with 
any  oxidizing  sulfide  ore  bodies  or  graphite  lenses. 

We  shall  use  our  model  of  a  vertically  layered  earth  to 
show  how  these  spontaneous  potentials  can  be  associated 
with  the  presence  of  a  contact  or  a  dike  through  the 
phenomenon  of  electrochemical  coupling.  It  will  be 
sufficient  to  use  a  two  dimensional  model  with  the 
z-axis  pointing  vertically  downward  along  one  of  the 
boundary  planes  and  the  x-axis  pointing  to  the  right 
(Fig.  3.8)  along  the  surface  of  the  ground.  We  shall 
treat  separately  the  vertical  contact  and  the  dike  pro¬ 
blems.  : 

Vertical  Contact  Problem  -  For  the  case  of  two  quarter 
spaces  separated  by  a  vertical  boundary  at  x«o  ,  we 
use  rectangular  coordinates  (  )  and  assume  that  the 

driving  potential  'Y  is  independent  of  y-coordinate . 

We  can  then  write  (see  Appendix  A) 


^  ^  R,e  Co* 

o 


X>  o 


X<o 


3.20 


Cn  A  £ 


3.21 
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Using  the  boundary  conditions,  we  get 


ft  -  io+ocvo  y,jX)  3.22 

P>1  »  -  3.23 


where  Rj.|  has  the  same  meaning  as  before  and  XiA)  is 
the  Fourier  Cosine  transform  of  the  driving  potential 
evaluated  at.  a:  =.  o  , 


Therefore,  for  potentials 


f  and 


we  get 


3.24 


pO  » 

f,  -  Co $Aa47l 


3.25 


3.26 


Equations  3.24,  3.25*  and  3*26  illustrate  once  more 
that  the  potential  -f  depends  only  on  the  value  of 
y(x/'%)at  the  boundary.'  In  this  case,  the  horizontal 
variations  af  'yr  are  unimportant  to  the  generation  of 
currents,  but  it  is  not  so  for  the  potential  which 

is  given  by 

f  foo  -  c  yc*,*) 


3.27 
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Furthermore,  the  above  equations  show  that  -f-  is  discon¬ 
tinuous  across  the  boundary,  while  ^  is  continuous. 

To  prove  that  is  continuous,  we  have  only  to  calcu¬ 
late  and  ^  on  both  sides  of  the  boundary: 

/ 

f/<V)  =  f ,<v)-  c.vm*  ifwi f.ofc-c.>y-c,/ 


oic^tr,  c, 

n+o*. 


XO/2) 


=  fAO  -  ^  ^ )  Cc,-cz)  y  _  cvy 

-  -  qrCt+<J7C^  yc<>,2) 


ffj  ■+  cq 


We  should  not  be  disturbed  by  the  discontinuity  of  f 
because  only  X  has  any  Physical  significance  and  can 
be  measured  (see  Appendix  B) . 

Again  if  the  boundary  does  not  cut  across  equi- 
potential  surfaces,  we  have 

Y(oj  %  )  =■  (a  constant) 

and  therefore 

■f|(*/2)  55  "2  0*  j  ")C 

=  -  -z0-  - 

Consequently  the  potential  4"  is  constant  in  each 
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medium  and  there  is  no  current  generated  as  predicted 
earlier. 

Application  to  Electrochemical  Diffusion  Potential  -  In 
order  to  get  some  numerical  results,  we  shall  apply 
equations  3.25  and  3.26  to  the  study  of  currents  generated 
across  a  contact  through  electrochemical  coupling.  This 
is  the  problem  commonly  encountered  in  drill  hole  elec¬ 
trical  logging.  Here  we  shall  assume  that  the  chemical 
potential  is  a  function  of  depth  Z  (it  may  also  be  a 
function  of  x-coordinate,  but  we  are  only  interested  in 
its  value  at  x  =  °  ) .  We  assume  that  its  z-dependence 
can  be  written 

-  kZ 

-  Ca/O  e,  3.28 

This  choice  of  u.  ,  although  arbitrary  is  applicable  to 
a  physical  problem  that  will  be  discussed  in  Chapter  V. 

The  chemical  potential  decreases  from  a  value  Jj.  = 
at  i  =.  o  to  /u=(iy4Q)el  at  a  depth  Z-1/^  which  might 
be  called  the  "effective"  depth  of  weathering,  after 
which  the  chemical  potential  remains  roughly  the  same. 

We  have  implicitly  assumed  in  3.28,'  that  at  great  depth 
the  chemical  potential  is  zero,  because  adding  a  constant 
value  does  not  affect  the  current  distribution  or  the 
potential  difference. 

Using  equations  3-24  and  3.25  we  get  (Erdelyi,  1954) 


/to)  = 2  .  .i 

^  ’T  V.K* 


3.29 


=  ka^„  J  -g~— c 3: 


30 


: a vk* 


Distance  (kx) 

FIGURE3.il  ELECTROCHEMICAL  POTENTIAL  <f>  (x)/  [(c,  -c2)(  A/i)]  FOR  A  VERTICAL 
CONTACT 
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This  solution  can  be  evaluated  in  terms  of  exponential 
integrals  with  complex  arguments.  For  i  =  o  ,  the  solu¬ 
tion  is  very  much  simplified  and  we  get  (Erdelyi,  195^0 


where 


3.31 


ciu*} 


and 


Mu)  --  - 


U. 


Figure  3.11  shows  the  behaviour  of  on  the  surface. 

It  is  Interesting  .to  note  that  far  from  the  common 
boundary  the  electrical  potential  ^  is  not  zero.  The 
two  media  are  at  different  potentials  and  the  maximum 
potential  difference  is 

=  Cc-Ci)  b/J-o  3.32 


For  the  current  flow  across  the  boundary  separating 
the  two  media,  we  have 


=  ^(we,occrQC^)<:^)  5  ^7 

J1x(o,2)  0+  ^»-X?i-Cv)(Kcr.)(^,)  [g  Et(w)_e  **£*(«)] 


3.33 
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where 

=  aa 

and 

<*o 

£*(«■)  -  -  J  ^ 


Vertical  Dike  -  A  more  realistic  case  is  that  of  a 
vertical  dike  of  conductivity  and  coupling 

coefficient  embedded  in  a  homogeneous  medium  (see 

Fig-.  3*8).  The  boundary  conditions  are 


^4,  _  -b  W 

~  -Qi- 


a4 


°7 


4" .  *■ 

iX,  1 

fe-ox 


11 


di  Z  =  Q 

. 


a.^ii  -  (r  lil 
•W  * 


V  ^z  =  ^C|'Cl)  X» 


3.34a 


3.34b 


3.34c 


The  two  media  on  both  sides  of  the  dike  are  identical  and 
therefore  =.  .  The  potentials  and  $ z  are 

of  the  following  forms 

=  J  ft,*  Cco  X*  <d*  3.35a 

f,  =  Ia  V**)^*^ 


3.35b 
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Substituting  from  3.3^b  and  3.34c  into  3.35a,  and  3.35b 
we  get  (see  Appendix  A) 


ft,  =  Ku**‘}(c.-cO 

V  -  ^0-Km)Cc<-cx) 


l-e  ^ 

vj» 

-U 

e 

■  Xp) 

1  0  -  Ad 

1- 

1 

%.(» 

and  finally  we  have 


0 O  _  ^ 

*. 3  i^K«X<r^)j  (-jf^cr)  j£oo£Axa>s>*d}  Oo  3>36a 

■ft®  ■jCl-l?*,-)CcrCj) j  (— — — J X&)  QI7^  *  <=0  -J  <?  x  <?*  3.36b 

Assuming  once  more  that  the  driving  potential  is  of  the 
form  3.28,  we  get  for  the  potentials  on  the  surface 


3.37a 


Wv)='~" ^ 5  fe-e*)W)c^)  j  H 


3.37b 
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Potential  is  plotted  in  Fig.  3.12  for  several  values 

of  and  f?t/  . 

Current  flow  across  the  boundary  is  given  by 


^^CW^)Cc-Cx)(K<r)C^)  j  - 


( \-c 


0 


Ut 


3.38 


and  is.  plotted  in  Fig.  3.13*  Near  the  surface  the  current 
flows  from  medium  II  into  medium  I  (for  cK),ct  )  and  the 
current  density  is  high.  At  a  depth  nearly  equal  to  the 
effective  depth  of  weathering  */«  ,  the  current  reverses 

its  direction  and  remains  that  way  for  increasing  depth 
but  the  current  density  goes  rapidly  to  zero. 

These  graphs  of  ^  and  for  our  hypothetical 
model  show  striking  resemblance  with  self -potential  field 
data  and  have  all  the  characteristics  usually  associated 
with  spontaneous  polarization  and  current  flow  around 
oxidizing  sulfide  veins. 

Section  IV  -  Random  Distribution  of  Inhomogeneities 

W$  shall  assume,  as  in  Section  II,  that  the  half- 
space  has  a  uniform  conductivity  c T  .  Therefore,  the 
potential  -f  is  due  to  distributions  of  dipoles  of 
strength  (c.-f,.)/  at  each  boundary.  Assuming  for  the 
present  that  the  coupling  coefficient  C  is  a  contin¬ 
uous  function  of  coordinates,  we  get 

=2  jj  ^ 

i  K 

+  2  jj  [c(%\  -C6d]  “  Ydc)**^* 

•  4  K  3.4! 

+  2  Si  [Cfifl-Cft)] 


0.6 


FIGURE  3.13  CURRENT  FLOW  Jlx(z)  /  [(c, -c2)(k<r,)(A^)J  ACROSS 
THE  DIKE  BOUNDARY 
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Clearly  we  can  write 

CCV)  -  c  (2-} 

and  similar  relations  for 
As  the  number  of  boundaries  tends  to  infinity,  the 
summation  over  index  Ji  must  be  replaced  by  an  integra¬ 
tion  and  we  finally  get 

2*4  =  jjj  [  VC.v(^)]y(2)dV  3.U3 

Equation' 3 -^3  is  the  exact  form  for  finding  -f  when  C 
is  known  as  a  function  of  coordinates.  Its  main  merit 
is  that  it  expresses  f  in  terms  of  a  volume  distribution 
of  dipoles . 

As  a  simple  model  of  random  medium,  we  shall  assume 
that  £> '  is  a  random  function  and  its  variation  depends 
only  on  the  X  -coordinate.  This  is  a  rather  drastic 
simplification,  but  it  can  be  used  as  a  model  for  a  series 
of  vertical  dikes  of  arbitrary  thicknesses  and  coupling 
coefficients.  Therefore,  we  can  write 

( yc)x  =  *i(x) 

We  shall  assume  that  0^  is  a  member  of  an  ensemble  for 
which  the  ergodic  theorem  holds,  i.e.  the  average  over  the 
ensemble! is  equal  to  the  average  over  the  x-coordlnate . 
Furthermore,  since  o(;(x)  is  the  gradient  of  the  coupling 
coefficient  C  (x)  ,  we  can  assume 


Ax- 


X- 


3.42 


^  and  i  variations  of  C 


(  *  ) 


o 


3.44 
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Substituting  3*^  in  we  get 

Ti  ^yfeodv'm  o 

This  means  that  the  induced  potential  oscillates  about 
a!  zero  average  value.  Taking  the  second  moment  to  find 
out  the  RMS  value  of  the  fluctuations,  we  get 


/  -x 


Rui 


3M 


In  the '.above  equation  o(;  represents  the  auto¬ 

correlation  function  of  and  we  shall  write  it 


«7  =  3.46 

i  ' 

Using  the  result 

f  <?JV  _  f  ^a1 _  _  ^ 

.  _-L  f?,J  “  J  [cxi^) vCb’-yjV (x-x)S  2-,z 

we  can  integrate  3.-^5  with  respect  to  Y  and  to  get 


_  1  (fff  ^  ( i^Lvk/)  U'_*\  Y(i'}  (*-x)YW  .  ,,  ,, , 


3.^7 


To  get  any  farther  we  need  to  make  specific  assump¬ 
tions  on  y(%)  and  ~\l^(lX-X"l)  .  Since  we  are  going  to 
use  the  model  to  predict  background  noise  due  to  electro- 
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chemical  diffusion,  we  shall  follow  Section  III  and  put 

Y<*)  =  e"** 


Making  a  shift  of  variables  and  xix  - — 

and  using  the  result  (Erdelyi,  1954) 

r*5  r 

fCH)  =  M  7T — Tdz  =  |sui(k^CrCt§)- 

O  i  J 

we  get'  !  , 

Si  ^OiVo  Kt?)r(tr) 


3.48 

Function  FCfc§,)  is  a  tabulated  function.  We  already 
have  used  this  function  in  Section  II  of  this  chapter. 

The  above  result  can  be  put  in  a  still  simpler  form  by 

making  the  change  of  variable  £  1 _ >  't+S"  •  We  then 

get 


^  rr2  k 
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where 


K  j~ (c|"J  F 


%p(^>)  = 


3.50 
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The  function  (k-irl)  is  the  autocorrelation  function 

of  F(X|)  which  represents,  except  for  a  constant 
coefficient,  the  potential  due  to  a  single  vertical  contact. 
According  to  3 .^9j.  the  rms  value  of  the  background  noise 
is  a  weighted  function  of  ^  (k-lri)  ,  the  weighting 

function  being  (irl)  .  The  autocorrelation  function 

L Tff  (klTl)  is  plotted  in  Fig.  3  •  1^  • 

To  get  specific  answers  for  ^  ,  one- has  to 

compute  from  the  field  measurements  of  wfx)  =  ^ 

No  such  extensive  measurements  have  been  undertaken  as  yet. 
Here,  we  shall  assume  that  the  autocorrelation  function 
li-  (in)  can  be  written 

*  otU 

4^  (in)  = 

This  choice  of  T^On  )  .>  although  ar 
justified  in  Chapter  V.  (AC)2"  is 
the  coupling  difference  AC*  between 
"t  is  the  average  thickness  of  the  d 

With  this  choice  of  l ^  (irl)  ,  w 

=  r  %,(<>) -<&?(/<*) 

TT1  <■  ki 

which  clearly  indicates  that  the  noise  level  is  a  function 
of  the  average  dike  thickness  t  as  well  as  the  chemical 
potential  /yx  and  average  variation  in  coupling  (AC) 
Application  of  these  equations  to  S.  P.  noise  will  be 
taken  ^in  Chapter  V. 


3.51 


bitrary,  will  be 
the  rms  value  of 
two  adjacent  dikes; 
ikes;  .  . 

e  get 


] 


3.52 
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The  function  (t(r/)  is  the  autocorrelation  function 

of  F(<$)  which  represents,  except  for  a  constant 
coefficient,  the  potential  due  to  a  single  vertical  contact. 
According  to  3-^9i.the  rms  value  of  the  background  noise 
is  a  weighted  function  of  (a/t/J  ,  the  weighting 

function  being  ^  (jrl)  .  The  autocorrelation  function 
(  W)  is  plotted  in  Fig.  3.1^. 

To  get  specific  answers  for  ,  one- has  to 

compute  ^Jin)  from  the  field  measurements  of  c<U)  =  ^ 

No  such  extensive  measurements  have  been  undertaken  as  yet. 
Here,  we  shall  assume  that  the  autocorrelation  function 
d'  (if'J  can  be  written 

^(m)  =  3,51 

This  choice  of  ,  although  arbitrary,  will  be 

justified  in  Chapter  V.  (AC)2-  is  the  rms  value  of 
the  coupling  difference  AC*  between  two  adjacent  dikes; 

■fc  is  the  average  thickness  of  the  dikes;  .  . 

With  this  choice  of  (tri)  ,  we  get 

_  LCtfCW-  F  %?(<>) -(PffMI  3.52 

jp  I  Jkt  J 

which  clearly  indicates  that  the  noise  level  is  a  function 
of  the  average  dike  thickness  t  as  well  as  the  chemical 
potential  /yu.  and  average  variation  in  coupling  (AC) 
Application  of  these  equations  to  S.  P.  noise  will  be 
taken  ^in  Chapter  V. 


CHAPTER  IV 


APPLICATION  TO  THE  DETECTION  OP  UNDERGROUND 
NUCLEAR  BLASTS 


Purpose 

Having  derived  the  equation  governing  the  distri¬ 
bution  of  potentials  induced  by  the  coupling,  of  flows 
and  studied  the  effect  of  geometry  on  these  coupling 
currents,  we  now  proceed  to  the  applications  of  the 
theory.  In  this  chapter,  the  potentials  generated  by 
the  electrokinetic  and  thermoelectric  couplings  are 
computed  to  show  the  applicability  of  these  couplings 
to  the  detection  of  an  underground  nuclear  blast  by 
its  electrical  aftereffects.  In  the  next  chapter,  the 
problem  of  spontaneous  polarization  near  inhomogeneities 
is  discussed  in  detail.  No  attempt  will  be  made  to 
cover  every  geophysical  application. 

Temperature  and  Pressure  Changes 

When  a  nuclear  charge  Is  detonated  underground,  some 
of  the  released  energy  serves  to  fracture  rocks,  produce 
seismic  waves,  etc.,  but  most  of  the  energy  remains 
confined  In  the  near  vicinity  of  the  blast  in  form  of 
heat  and  serves  to  raise  the  temperature  of  the  pore 
fluids  and  create  steam  under  pressure  in  equilibrium 
with  the  fluids.  The  earth  acts  essentially  as  an  infin¬ 
itely  large  fluid  reservoir,  and  therefore  the  steam 
temperature  never  rises  above  the  boiling  point  temperature. 
This  point  is  confirmed  by  drill  hole  temperature  measure¬ 
ments  made  a  few  days  after  the  Rainier  shot  (Johnson  and 
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8S°C  in  the  blast  area.  Incidentally,  this  is  nearly 
equal  to  the  boiling  point  temperature  of  water  for  the 
atmospheric  pressure  at  those  devations  (The  mesa  in 
which  the  charge  was  detonated  is  several  thousand  feet 
above  the  sea  level) . 

If  the  explosion  is  totally  contained  and  no 
extensive  fracturing  takes  place,  the  only  way  the  steam 
can  escape  is  by  forcing  the  surrounding  fluid  out. 

Under  these  conditions,  the  high  pressure  source  should 
exist  for  a  long  time  as  shown  in  Chapter  I.  But  it  is 
expected  that,  even  though  the  blast  may  take  place  at 
depths  of  several-  thousand  feet,  extensive  fracturing 
should  take  place  which  would  increase  the  permeability 
of  the, medium,  and  speed  up  the  release  of  pressure. 

Also  if  the  gas  should  channel  its  way  out,  we  could 
expect  an  even  faster  release'  of  pressure.  The  actual 
situation  is  difficult  to  predict. 

When  the  shot  is  not  buried  too  deeply,  the  fractur¬ 
ing  may  open  up  fissures  reaching  to  the  surface,  and 
thus  release  the  excess  steam  until  the  pressure  inside 
the  blast  cavity  is  reduced  to  the  atmospheric  pressure. 
But,  as  we  have  already  pointed  out  in  Chapter  I,  even 
in  this  case  the  pressure  effect  would  not  be  lost, 
because  the  excess  of  the  hydrostatic  pressure  over  the 
cavity  pressure  would  cause  a  backflow  of  fluids  into  the 
cavity.  For  a  shot  buried  at  200  meters,  this  excess 
back  pressure  can  be  of  the  order  of  20  to  '50  atmospheres 

While  the  temperature  picture  we  presented  here  has 
been  confirmed  by  drill  hole  temperature  measurements 
after  the  Rainier  shot,  no  direct  evidence  exists  to 
confirm  our  pressure  model.  But  there  is  indirect 


evidence  to  suggest  that  our  pressure  reversal  model  can 
be  applied  In  certain  cases.  ~  ^  ^  - 

MSffiffflfPoverTa r dvark  shot  showed  an  abrupt  change  in 
sign  that  could  be  easily  explained  by  our  model. 

Chemical  reactions  that  follow  the  blast  could  possibly 
add  to  the  complexity  of  the  picture.  They  are  little 
understood  except  that  the  time  scales  involved  are  very 
long,  and  therefore,  the  role  of  electrochemical  diffusion 
as  a  source  of  currents  can  not  be  evaluated  at  this 
point.  However,  the  temperature  and  pressure  effects  are 
understood  enough  to  attempt  to  predict  the  magnitudes  of 
thermoelectric  and  electrokinetic  potentials  to  be  expected, 
and  compare  them  with  the  potentials  measured  at  some  of 
the  underground  nuclear  detonations . 


Thermoelectric  Effect  of  a  Bomb 

As  we  have  already  shown,  the  temperature  in  the 
blast  area  can  not  be  above  the  boiling  point  temperature 
of  water  at  the  prevailing  pressure.  Let  us  assume  that 
the  pressure  is  atmospheric  and  therefore  this  boiling 
point  temperature  is  about  100°C.  Assuming  an  ambient 
temperature  of  20°C,  we  get  a  temperature  difference 
AT  3t  80 *C  .  If  the  shot  were  completely  contained,  the 
pressure  would  build  up,  and  consequently,  this  estimate 
of  80°C'  would  be  far  too  low. 

The  thermoelectric  coefficients  of. rocks  (or  rather 
the  difference  between  the  coefficients  for  rocks  and 
that  for  water)  are  seen  from  Table  2.4  to  have  any  value 
between  0  and  I.36  mv/°C,  with' typical  values  of  .2  to  .5 
mv/°C.  Therefore,  a  typical  value  for  the  difference  of 
coupling  (  Cy-Ct,  )  would  be  about  .2 — .3  mv/’C.  This  is 
only  an  order  of  magnitude  and  it  could  be  as  large  as 


one 


millivolt  per  degree  or  as  little  as  zero. 

To  get  the  geometric  factor  from  Fig.  3*4,  let  us 
assume  that  a  20KT  bomb  is  exploded  inside  a  half  space 
at  a  depth  of  100  feet  below  a  homogeneous  layer  as  shown 
in  Fig.  3.2.  If  uhe  thickness  of  the  overlying  layer  is 
300  feet  and  the  rau.  us  of  the  high  temperature  cavity 
is  about  l80  feet,  the  radius  0.  of  the  high  temperature 
zone  intercepted  at  the  "boundary  between  the  layer  and  the 
half  space  would  be 

O-  =  \[l  80t-|oolt-  3s  ISO  f+ 

and  the  ratio  (^)=  •  5 •  Assuming  a  conductivity  contrast 

,  and  using  the  data  of  Fig.  3 .4,  equation  3*17 

gives  us 


Af  ^  =•  (o. 3  m>//^d)(©0*c)(.|l)  =  2.C,  m\j  4.1 

The  factor  (.11)  in  4 . 1  is  the  geometric  factor  and  it 
indicates  that  although  "the  induced  thermoelectric 
potential  is  24  milllvol-ts,  the  maximum  potential  observed 
at  the  surface  is  only  about  eleven  per  cent  of  that  value. 

Thi3  estimate  of  2.6  mv  is  too  small  to  be  easily 
detectable.  A  refinement  of  the  present  exploration 
techniques  along  with  a  thorough  study  of  the  nature  of 
the  background  noise  would  be  necessary  to  develop  tech¬ 
niques  for  detecting  such  potentials. 

This  small  voltage  indicates  that  the  thermoelectric 
coupling  is  not,  in  general,  very  important  in  producing 
electric  potentials  from  a  nuclear  blast.  Under  very 
favorable  conditions,  i.e.,  half  space  much  more  conductive 
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uicui  ...... ,  wW; , .  -ii^vryv -A-fi t*.  cuc-e-  vi'-  'couplings 

(tf-Ct)  >  O.ftHv^rJ  a  ratio  (t)-^0  >  and  a  temperature 
difference  kT^Zoo'C  i  we  would  get 

£fmay  -  (o^^/6c)(-2o°‘c)f.s)  mu  2^2 

which  is  a  quite  large  potential  difference  and  can  be 
easily  measured  with  the  present  techniques  and  instruments 
However,  coupling  differences  of  osn  ^°c  and  temperature 
.differences  of  200 °C  can  be  considered  unreasonably  high 
in  most  cases,  and  therefore  thermoelectric  effects  of 
50  mv  Should  not  be  expected. 

We  conclude  this  section  with  a  comment  on  Fig.  3.4 
showing  the  geometric  factor  for  the  thermoelectric  effect. 
Over  the  shot,  the  geometric  factor  f/r)/ft.V<:o'7;’] is  greater 
for  positive  values  of  K?.\  (i.e.,  half  space  more  con¬ 

ductive  than  the  layer);  but,  at  large  distances,  the 
reverse  is  true.  This  means  that — assuming  thermoelectric 
potentials,  are  large  enough  to  be  detected — a  bomb  can 
be  more  easily  detected  if  the  overlying  layer  is  a  rela¬ 
tively  good  conductor,  because  the  effect  persists  over 
a  larger  area . 

Electrokinetic  Effect  of  a  Bomb 

Streaming  potential  measurements  shown  in  Table  2.3 
indicate  that  there  are  large  variations  in  the  electro- 
kinetic  coupling  coefficients  of  the  rocks,  with  the 
smaller  coefficients  corresponding  usually  to  more  tightly 
packed  rocks.  Table  4.1  summarizes  these  results: 
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. **E m&r  -A..J  ^  ^%miaTy~oIr>S,t?eaming/Potential  Measure¬ 


ments  Made  on  Rock  Samples 


Number 

Rock  Type  of 

Samples 


Sedimentary  Rock, 

Colorado  Plateau  9 

Dakota  Sandstone  6 

Altered  Volcanics, 

Arizona  10 

Metamorphlc  Rocks, 

Ontario  5 

Ophytic  Trap  Recks, 
Michigan  2 

Quartz-Latite  Porphyry, 

New  Mexico  4 


Streaming  Potential  (mv/atm.) 


Average 

Value 

Standard 

Deviation 

13.8 

13.5 

7.3 

7-5 

8.4 

7-5 

8.9 

5.2 

4.3 

0.4 

.40 

.16 

We  already  have  discussed  briefly  what  happens  when 
a  large  charge  is  detonated  underground.  The  pressure 
bulld-Up  that  follows  is  influenced  by  several  factors,  and 
the  resulting  pressure  picture  is  difficult  to  predict 
theoretically. 

If  the  blast  were  completely  contained  and  no  gases 
were  allowed  to  escape,  the  pressure  Inside  the  cavity 
would  be  definitely  greater  than  the  ambient  pressure.  A 
simple  calculation  will  give  us  the  order  of  magnitude  of 
this  excess  pressure.  Assume  that  the  rocks  are  saturated 
with  water,  and  the  volume  increase  due  to  outflow  is 
negligible.  Therefore,  any  steam  produced  will  have  the 

O 

same  specific  volume  as  water  (18  cirr/mole).  Using  a 
temperature  of  373°K  and  assuming  that  steam  at  this  temper¬ 
ature  behaves  like  an  ideal  gas,  we  get 


p  -  51  -  g-*2-*3"75 _  £  111  XID4  ne«)WAl 

^  Ift  X  to-6 

or  ±  -  1700  atmospheres . 

This  estimate  will  be  too  high,  partly  because  rock 
porosity  is  always  greater  for  gases  than  for  fluids  and 
the  resulting  steam  occupies  a  larger  volume  than  the 
original  fluid,  and  partly  because  steam  at  100°C  is  far 
from  being  an  ideal  gas,  and  its  pressure  will  be  far 
below  that  predicted  by  the  perfect  gas  law. 

For  a  particularly  simple  model  of  a  blast  creating 
extensive  cracks  and  allowing  the  gases  to  escape  into  the 
air,  the  pressure  inside  the  cavity  would  be  the  atmos¬ 
pheric  pressure.  If  the  center  of  the  blast  is  located 
at  a  depth  of  400  feet  below  the  surface,  the  back  flow 
pressure  due  to  gravity  would  be  about  12  atmospheres. 

Let  us  compute  the  streaming  potential  for  back 
flow  conditions,  using  a  value  of  10  mv/atm.  for  the 
difference  in  coupling  (  C,-~  )„  We  assume  again  that 

the  blast  area  has  a  radius  a  =  180  ft.,  and  is  buried 
at  a  depth  ^  =  400  feet;  the  overlying  layer  has  a 
thickness  =  300  ft.,  and  the  conductivity  contrast 
is  zero.  We  then  have  *6,  (—)  =  and  from 

Fig.  3.5*  we  get  an  approximate  value  of  0.5  for  the 
geometric  factor  ^(r) ]  right  over  the  shot 
(  rso  ).  Thus  the  maximum  potential  difference  is 


kfma*  =  0°  ™ - ) (°-0 ( i 2- J (a.  s')  =  34  ™\/  4.3 
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This  is  a  sizeable  anomaly  and  could  be  detected. 

Under  very  favorable  conditions,  the  backflow 
pressure  can  be  much  higher;  pressures  of  50  atm.  or 
better  are  possible.  Using  a  drastically  large  value 
of  20  mv/atm.  for  the  coupling  difference,  we  get 

(2  O  mM Mr*  )  (oi)  (So  oJcr*. )  (o s')  ■s  3oo  rAv 

4.4 

This  represents  a  very  large  anomaly  and  should  be  consid¬ 
ered  as  an  upper  limit  of  the  potentials  to  be  expected. 

The  actual  potentials  measured  at  some  of  the  underground 
nuclear  detonations  are  only  of  the  order  of  a  few 
millivolts  (See  Pig.  4.2,  4.3,  and  Table  4.2)., 

P.rom  these  results  and  those  of  the  previous  section 
we  should  not  conc3.ude  that  the  streaming  potential  is 
necessarily  the  more  important  effect.  That  would  depend 
on  the  local  geologic  conditions  which  could,  for  instance, 
produce  large  thermoelectric  and  small  electrokinetic 
coupling  differences  at  the  same  time.  It  is  important 
to  note  that  only  differences  in  coupling  are  involved. 

It  is  possible  to  have  two  layers  with  large  coupling 
coefficients  and  yet  observe  no  streaming  potential  effect 
because  the  difference  in  coupling  is  small.  Another 
possibility  is  to  have  a  situation  in  which  a  thin  layer 
separates  two  media  of  comparable  coupling  effects.  In 
this  case  the  differences  (  c(_c,  )  and  (  )  are  of 

opposite  sign  and  therefore  the  induced  potentials  tend 
to  cancel  each  other  (see  equation  3-19).  The  above 
results  indicate,  however,  that  in  principle,  the  electro- 
kinetic  coupling  can  be  the  more  important  of  the  two. 

At  least,  the  results  of  Nevada  underground  tests  tend  to 
confirm  this  view. 


S*IOA!H!W 


FIGURE  4.!  ELECTRIC  POTENTIAL  DECAY  FOR  AARDVARK  SHOT,  PROJECT  NOUGAT 
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Nevada  Underground  Detonations 

As  part  of.  an  investigation  of  the  feasibility  of 
using  surface  electrical  measurements  to  detect  under¬ 
ground  nuclear  detonations,  extensive  field  measurements 
were  undertaken  in  conjunction  with  Project  Gnome  and 
some  events  of  Project  Nougat  (Glaser  and  Kruszyna,  1962; 
Kruszyna,  1962).  There  were  noticeable  potential  changes 
of  10  to  30  mv.  in  many  cases.  The  results  from  these 
shots  are  summarized  in  Table  4.2. 

The  observed  potentials  at  the  Gnome  event  showed 
a  sharp  rise  over  a  period  of  36  hours  followed  by  a 
gradual  decay  toward  the  ^re-shot  noise  level.  The 
decay ^was  complete  after  eight  days.  If  our  pressure 
model ^is  applicable  here,  it  would  suggest  the  no  degassing 
took  place  and  the  observed  effects  were  due  to  streaming 
potentials  created  by  out-flowing  fluid. 

The  results  of  the  Hardhat  event  show  as  sharp  poten¬ 
tial  rise  over  a  period  of  two  days  followed  by  an  abrupt 
drop  within  a  few  hours,  after  which  the  potential  started 
to  rise  again  above  the  pre-shot  level,  and  it  was  still 
rising  seven  days  after  the  detonation  when  the  recording 
ceased,.  The  potential  reversal  in  this  case  could  be 
explained  by  our  model  of  pressure  release  but  it  would 
be  difficult  to  explain  why  the  potential  continued  to 
rise  again  above  the  ambient  level. 

The  Ardvark  event  was  the  most  conclusive  of  the  six 
tests.  A  large  and  sharp  potential  rise  of  about  24  mv. 
was  followed  by  an  abrupt  reversal  to  a  value  well  below 
the  ambient  level  (Pig.  4.1).  after  which  the  potential 
decayed  over  a  period  of  45  days  to  its  pre-shot  level. 

This  behaviour  fits  very  well  our  model  of  flow  reversal 
caused  by  degassing.  The  A.  C.  noise  level  was  relatively 


4 
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Table  4.2  -  Summary  of  Potential  Measurements  from 

Projects  Gnome  and  Nougat 


Event' 

Gnome 

Hardhat 

Black 

Aardvark 

Eel 

Marshmallow 

Yield 

Low 

5  KT 

Low 

5  KT 

Low 

Medium 

37  KT 

Low 

Low 

Depth 

(feet) 

1200' 

950' 

1434' 

1200' 

Magnitude  of 
Potential  use 
(mv) 

15 

40 

10 

24 

3? 

Magnitude  of 
Potential 
reversal  (mv) 

30 

30? 

35 

10 

i3? 

Decayed  to  Pre¬ 
shot  Value? 

yes 

no 

9 

yes 

yes 

9 

Approximate 

Decay  time 
(days) 

8 

? 

50? 

45 

16  ? 

? 

Noise  Level 

Low 

Moder¬ 

ate 

High 

Low 

Low 

Moder¬ 

ate 
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low  during  this  event  and  good  enough  records  are  available 
to  attempt  a  quantitative  interpretation  of  the  results. 

It  is  possible  that  the  good  quality  of  the  records  was 
due  to  the  large  size  of  the  blast.  This  was  the  only 
medium  yield  shot  of  the  series. 

In  the  case  of  the  Black  event,  even  though  relativelt 
large  potential  changes  were  observed,  no  definite  conclu¬ 
sion  can  be  drawn  due  to  the  large  variations  in  the  noise 
level. 

The  Eel  event  did  not  show  any  potential  rise,  but 
rather  the  potential  dipped  some  10  mv.  below  the  background 
level  and  gradually  recovered  in  about  two  weeks. 

In  the  marshmallow  event,  a  very  modest  potential  rise 
of  3  mv.  occurred  over  a  period  of  two  days.  A  potential 
reversal  followed,  and  at  the  end  of  3  days,  the  observed 
potential  was  some  13  mv.  below  the  pre-shot  level.  The 
recording  was  stopped  after  seven  days  making  it  impossible 
to  decide  whether  the  readings  were  due  to  blast  induced 
potential  changes  or  erratic  telluric  currents. 

Aardvark  Shot 

Aardvark  event  was  a  37  +  7  KT  detonation  buried  at 
a  depth  of  1434  ft.  inside  a  tuff  layer  overlain  by 
alluripm.  The  thickness  of  the  alluvial  layer  is  not  known 
exactly  but  it  is  of  the  order  of  1000  ft.  Five  potential 
measuring  electrodes  were  located  radially  at  250  ft. 
intervals,  the  nearest  being  located  about  500  ft.  from 
ground:  zero.  >  All  potential  measurements  were  made  rela¬ 
tive  tp  a  reference  electrode  R  located  at  4000  ft. 
from  ground  zero.  Measurements  were  made  starting  4  days 
beforej  and  continuing  to  '53  days  after  the  blast.  A 
potential  change  between  each  of  the  five  electrodes  and 
the  reference  electrode  was  observed.  A  sharp  potential 
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rise  was  followed  In  each  case  by  an  abrupt  drop  to  a 
value  below  the  natural  D.  C.  level.  Figure  4.1  shows  the 
potential  difference  between  electrode  1  (the  nearest  to 
the  shot)  and  the  reference  electrode  R.  .  The  line 
marked  "D.  C.  Level"  is  an  estimate  of  the  background 
potential  that  was  not  completely  bucked  out.  The  broken 
line  indicates  the  measured  potentials,,  while  the  heavy 
line  represents  the  probable  signal  due  to  the  blast. 

The  large  potential  reversal  and  the  decay  time  of  about 
45  days  are  two  very  important  features  that  will  be 
discussed  in  the  next  section. 

The  measurements  from  the  other  electrodes  have  the 
same  general  features,  although  the  magnitudes  of  the 
observed  potentials  are  systematically  smaller.  In  each 
case  there  was  a  positive  potential  peak  followed  abruptly 
by  a  negative  potential  peak  and  a  gradual  decay  over 
seven  :or  eight  weeks.  Figure  4.2  shows  the  positive  peaks 
by  a  series  of  bars.  The  length  of  each  bars  indicates 
the  uncertainty  in  the  interpretation  due  to  the  presence 
of  noise.  Fig.  4.3  presents  similar  results,  except  that 
here  we  have  plotted  the  observed  negative  peaks.  Both 
graphs; show  clearly  that  the  induced  potential  decreases 
with  increasing  distance  from  ground  zero  as  predicted  by 
the  theory.  However,  they  differ  drastically  among  them¬ 
selves;  jin  exhibiting  different  rates  of  decrease. 

Tljie  time  scale  of  45  days  and  the  sign  reversal  as 
shown  in  Fig.  4.1  suggest  that  electrokinetic  coupling  is 
responsible  for  the  observed  electric  potentials.  We 
showed  in  Chapter  I,  that  the  time  scale  for  thermal  effect 
is  always  of  the  order  of  years,  while  shorter  time  scales 
could  be  expected  for  the  pressure  effects.  A  time  scale 
of  45  days  Is  very  compatible  with  the  calculations  made 
in  Chapter  I  for  the  fluid  backflow  case,  but  it  is  com- 
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Figure  4.3*  Observed  Negotive  Potentiol  Peok  vs.  Distonce  from 
-  Ground  Zero  for  Aordvork  Shot 
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pletely  unrealistic  for  thermal  dissipation.  The  sign 
reversal  is  still  a  more  positive  indication  of  the  impor¬ 
tance  of  the  streaming  potential  in  this  case.  The 
temperature  effect  cannot  be  expected  to  reverse,  but  if 
che  effect  were  due  to  electrokinetic  coupling,  the  sign 
reversal  could  be  explained  by  a  reversal  in  the  direction 
of  flow  following  a  pressure  release  due  to  the  outgassing 
of  the  cavity.  Thus  in  Pig.  4.1,  the  positive  peak  would 
be  interpreted  as  being  due  to  streaming  potentials 
induced  by  the  steam  pressure  forcing  the  fluids  out.  But 
in  a  matter  of  hours,  all  of  the  excess  steam  would  have 
leaked  to  the  surface,  resulting  in  a  slower  back  flow 
which  ...would  be  responsible  for  the  negative  potentials. 

Quantitative  Interpretation  of  the  Aardvark  Results 

As  we.  have  already  seen,  the  general  features  of 
the  Aardvark  shot  can  be  explained  on  the  basis  of  a 
reasonable  model  involving  electrokinetic  coupling  and 
flow  reversal.  We  wish  to  examine  more  carefully  the 
check ^between  the  field  data  and  the  theoretical  stream¬ 
ing  potentials  for  our  simple  model  assuming  one  layer 
over  half  space  geometry  (Figures  3.2  and  3-5).  In  par¬ 
ticular,  we  would  like  to  predict  the  depth  of  the  shot  and 
the  thickness  of  the  overlying  layer. 

One  unusual  feature  of  the  observed  data  is  the  rate 
of  decrease  of  the  potentials  with  increasing  distance. 

The  maximum  positive  peaks,  shown  in  Fig.  4.2,  decrease 
rapidly  with  increasing  distance  indicating  a  shallow 
source,  while  the  negative  peaks,  shown  in  Fig.  4.3, 
decrease  much  more  slowly  indicating  a  more  deep  seated 
source.  Furthermore,  the  positive  potential  peaks  last 
for  only  a  few  hours  and  therefore  represent  transient 
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conditions  of  very  short  duration.  The  plot  of  the  nega¬ 
tive  peaks  would  represent  the  more  stable  conditions 
following  the  short  period  transient  behaviour,  and  thus, 
would  be  more  reliable  for  quantitative  interpretation. 

The  positive  potentials  could  be  caused  by  upward  channel¬ 
ing  of  the  steam.  The  blast  could  have  opened  up  many 
cracks,  allowing  the  steam  to  push  its  way  upward  and 
thus  produce  a  shallow  and  short-lasting  streaming  poten¬ 
tial  effect.  The  reversed  potentials  that  followed  would 
represent  the  streaming  potential  created  by  the  backflow 
of  fluids  into  the  blast  cavity. 

The  theoretical  streaming  potentials  for  one  layer 
over  half  space  (Fig.  3*5)  have  been  replotted  in  Fig.  4.4 
using ;'the  depth  of  the  shot  "ft  rather  than  the  thickness 
of  the  layer  d  as  the  constant  parameter.  These  plots 
show  a  great  resemblance  in  their  rates  of  fall-off  with 
horizontal  distance.  The  thickness  of  the  overlying 
layeri'affects  the  magnitude  of  the  potentials  slightly, 
but  it  does  not  affect  the  shape  of  the  curves.  This 
indicates  that  the  depth  estimates  will  not  be  greatly 
influenced  by  the-  particular  geometry  that  is  assumed, 
the  paramount  factor  being  the  depth  of  the  shot. 

The  observed  data  from  the  Aardvark  event  have  been 
superimposed  on  each  of  the  three  plots  in  Fig.  4.4.  The 
circled  points  represent  the  potential  estimates  read 
from  Fig.  4.3.  It  can  be  seen  that  field  data  and  the 
computed  curves  fit  as  well  as  the  accuracy  of  the  data 
permits.-  The  match  of  the  field  data  with  the  streaming 
potential  curve  for  gives  "ft  =  1400  ft.  and 

d  =  1400  ft.  Using  the  plot  for  ,  we  get 

-ft  =  1650  ft.  and  d  =  825  ft.;  and  finally,  the  plot 
for  i^)~  3  gives  ft  =  2000  ft.  and  d  =  670  ft.  The 
Aardvark  shot  was  actually  detonated  in  a  tuff  layer  at  a 
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Figure  4.4*  Match  of  Observed  Data  with  Theoreticoi  Streaming  Potential 
Curves  for  One  Layer  over  Half  Space 
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depth  of  1434  feet.  The  thickness  of  the  overlying 
alluvium  is  not  known  exactly,  but  it  is  of  the  order  of 
1000  feet.  Considering  the  accuracy  of  the  data,  these 
estimates  of  $  and  d  a]re  not  too  far  off.  One 
disturbing  fact  is  that  these  estimates  of  the  depth  tend 
to  be  too  high.  Estimates  that  fall  short  of  the  actual 
depth  can  be  due  to  the  rock  fracturing  and  upward  funnel- 
ing  that  follow  the  detonation  and  have  the  effect  of 
moving  up  the  center  of  the  blast  after-effects.  On  the 
other  hand,  estimates  that  are  too  high  could  hardly  be 
explained,  and  would  be  looked  upon  with  suspicion. 

Estimates  of  the  coupling  coefficients  suggest  that 
our  pressure  model  for  the  :Aardvark  shot  is  very  reasonable.. 
According  to  equation  3.18.-  the  streaming  potential  is 
given  :by 

f  -  h.-MCtor. 4 (r, 4-5 

where  !?«)  is  the  dimensionless  "geometric  factor" 

plotted  in  Fig.  3.5*  and  Tp  is  the  driving  pressure. 

The  geometric  factor  GCn^e,,)  is  a  positive  number, 
therefore,  the  sign  of  the  potential  f  is  given  by  the 
sign  of  the  product  (c,  -  .  in  our  model  of  the 

Aardvark  shot,  we  said  that  the. negative  potentials  were 
caused  by  a  backflow  pressure.  Therefore,  in  our  case, 
is  negative  and  the  coupling  difference  (q-c,)  is  positive, 
indicating  that  the  overburden  has  a  larger  coupling 
coefficient  than  the  tuff  layer.  Indeed,  it  is  found  that, 
in  general,  the  loosely  packed  rocks  inhibit  higher  coup¬ 
ling  effects.  Using  a  shot  depth  of  1400  ft.  (the  top 
curve  in  Fig.  4.4),  we  get  for  the  geometric  factor  C n-cuz.  , 
and  from  Fig.  4.1,  we  get  kf  »  tw  .  The  radius  ct  of  the 
area  affected  by  a  37  KT  blast  is  around  200  ft.  giving  a 
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ratio  (^  =  o-)+  .  Assuming  that  at  1400  ft.  of  depth, 

the  hydrostatic  pressure  is  about  42  atm.,  we  get  for  the 
coupling  difference 


caw 


; -  «  3. rr^v^vyi.  ^-6 

<S>hX42)(o.u) 


This  estimate  of  3.3  mv/atm.  is  very  realistic  and  provides 
an  additional  check  of  our  pressure  model. 


Summary  of  Chapter  IV 

While  a  host  of  irreversible  effects  may  follow  the 
detonation  of  an  underground  nuclear  blast,  theoretical 
estimates  of  the  expected  potentials  as  well  as  the 
observed  data  indicate  that  the  streaming  potential  is  the 
most  important  after  effect.  It  can  generate  potentials 
of  sufficient  magnitude  to  be  detected  by  present  explora¬ 
tion  techniques. 

Comparison  between  the  field  data  and  the  theoretical 
streaming  potentials  for  one  layer  over  a  half-space  indi¬ 
cate  the  feasibility  of  quantitative  interpretation  of  the 
field  measurements.  The  match  between  the  observed 
potentials  and  three  different  theoretical  curves  gave  a 
fair  estimate  of  the  depth  of  the  shot,  the  error  being 
attributed  mostly  to  the  large  noise  level.  Some  of  the 
negative  potential  peaks  read  from  the  field  records  might 
be  off  by  as  much  as  100  per  cent. 

A  simple  model  of  pressure  release  and  flow  reversal 
was  found  to  explain  all  of  the  features  of  the  potentials 
created  by  the  Aardvark  shot.  However,  this  model  is  not 
expected  to  be  universally  applicable.  Indeed,  the  data 
recorded  at  Gnome  shot  did  not  show  any  reversal  of 
potential.  More  field  measurements  taken  with  improved 


-  79  - 


recording  techniques  are  necessary  before  definite  conclusions 
can  be  drawn  on  the  applicability  of  the  streaming  potential 
model  to  on  site  detection  of  underground  explosions. 


CHAPTER  V 


APPLICATION  TO  SPONTANEOUS  POLARIZATION  CURRENTS 

Spontaneous  polarization  currents  or  self-potentials, 
refer  to  natural  earth  currents  that  are  relatively  stable 
in  time — both  in  direction  and  magnitude.  They  are 
commonly  associated  with  the  weathering  of  sulfide  ore 
bodies,  the  movement  of  underground  water,  variations  in 
rock  properties  across  geologic  contacts,  the  bioelectric 
activity  of  plant  roots,  and  a  host  of  lesser  important 
phenomena  such  as  corrosion  of  pipe  lines.  Electrochem¬ 
ical  action  arising  from  concentration  differences  of 
electrolytes  in  contact  with  metallic  conductors,  spon¬ 
taneous  polarization  of  two  dissimilar  metals  in  contact 
with  solutions,  electrofiltration,  and  topographic  effect 
due  topthe  normal  atmospheric  electric  gradient  have  been 
proposed  to  be  the  most  important  mechanisms  of  generation 
os  these  currents . 

The  electrochemical  and  electrokinetic  mechanisms  of 
self-potentials  have  been  discussed  in  detail  by  various 
authors  (Heiland,  19^6;  Mounce  and  Rust,  19^5;  Sato  and 
Mooneyy.  i960),  but  no  serious  attempts  have  been  made  to 
evaluate  the  effect  of  geometry  on  the  magnitude  of  the 
induced  potentials.  We  shall  attempt  to  apply  the 
results  of  Chapter  II  and  III  to  an  evaluation  of  the  role 
of  various  couplings  in  producing  self-potentials  as  well 
as  an  estimation  of  the  effect  of  geometry  on  these 
potentials.  In  particular,  we  shall  examine  the  importance 
of  the: geothermal  gradient  in  producing  large  scab  S.  P. 
currents  through  thermoelectric  coupling,  the  role  of 
electrochemical  coupling  in  producing  small  scale  S.  P. 
currents,  and  the  effect  of  geometry  on  the  magnitude  of 
S.  P.  noise  produced  by  electrochemical  coupling. 
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Section  I  -  S.  P.  Currents  Caused  by  Thermoelectric  Coupling 

In  Section  III  of  Chapter  III,  we  proved  that  when 
vertical  boundaries  cut  across  the  equipotential  surfaces 
associated  with  a-  depth-dependent  driving  potential  , 
spontaneous  polarization  currents  are  generated.  The 
general  features  of  these  currents  are  shown  in  Fig.  3.11, 
3.12,  and  3*13.  It  is  believed  that  the  same  type  of 
currents  could  be  Induced  in  the  earth's  crust  where  large 
lateral  changes  of  rock  properties  and  a  vertical  tempera¬ 
ture  gradient  exist  at  the  same  time.  Indeed,  it  has  been 
proposed  in  the  past  that  , thermoelectric  currents  through 
the  mantle  and  the  crust  might  be  responsible,  at  least 
in  parfi,  for  the  earth's  magnetic  field. 

In  computing  the  thermoelectric  currents  in  the 
earth, '-.we  shall  limit  ourselves  to  currents  generated  in 
the  crust,  because,  in  the  mantle  the  conditions  are  vastly 
different.  The  experimental  data  on  thermoelectric 
coupling  coefficients'  listed  in  Table  2.4  are  for  porous 
rocks  where  the  conductivity  is  mainly  electrolytic,  and 
the  coupling  coefficients  for  the  rocks  in  the  mantle, 
where  the  mode  of  condition  iS  very  different,  are  not 
known.  ;• 

It  is  well  known  that  there  is  a  geothermal  gradient 
of  about  30°C/km  which  can  be  assumed  to  exist  at  least 
as  far Mas  the  base  of  the  crust.  This  will  give  an  appre¬ 
ciable!  i temperature  difference  of  some  1000°C  between  the 
top  and  the  bottom  of  the  crust  which  we  shall  assume  to 
be  about  35  km  in  thickness.  Such  a  large  temperature 
difference  can  create  sizeable  potential  difference  over 
distances  of  several  kilometers. 

Lqt  us  compute  the  magnitude  of  the  potential  differ¬ 
ence  along  the  surface,  assuming  that  the  lateral  crustal 
changes  in  coupling  properties  can  be  represented 
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approximately  by  a  vertical  discontinuity.  Furthermore, 
we  shall  assume  that  the  dependence  of  temperature  on 
depth  can  be  approximated  by  an  equation  of  the  form  3.28, 
where  now  the  "depth"  ^  is  of  the  order  of  20  to  30  km. 

We  can  then  use  the  equations  3-31  and  3.33  derived  for  a 
two-dimensional  vertical  contact  model.  The  conductivity 
contrast  Rj.i  can  be  put  to  zero  without  loss  of  generality. 
As  it  can  be  seen  from  Fig.  3.H>  the  potential  difference 
across  the  contact  is  very  little  affected  by  varia¬ 
tions  in  conductivity.  This  difference  of  potential  t 
for  two  points  on  opposite  sides  of  the  contact  and  a  dis¬ 
tance  (  Kx  )  =  1.0  away  from  it,  is  found  from  Fig.  3.11 
to  be  about  0.6  and  we  have  to  multiply  this  value  by 
the  coupling  difference  (  C\  -  C*.  )  and  the  temperature 
difference  AT  to  get  the;, value  of  the  potential  differ¬ 
ence  ini  millivolts .  Using  a  typical  value  of  0.2mv/°C 
for  the  coupling  difference,  we  get 

=  (o.z  mv/°c)  (\ooo°c)  (o.G)  =  \2 .0  mv  5.1 

This  rather  large  potential  difference  is  developed, 
as  we  said  above,  over  a  dimensionless  distance  (  Kx)  -  2.0; 
using  the  value  of  25  km  for  '/*  ,  we  get  for  the  actual 
distance  X  =  50  km.  Therefore,,  the  potential  drop  is  no 
more  than  2  or  3  mv/km.  From  the  point  of  view  of  the 
mining  geophysicist,  this  would  represent  a  rather  small 
noise  that  would  not  influence  his  studies. 

It  is  obvious  that  the  contribution  of  the  thermo¬ 
electric  currents  in  the  crust  to  the  earth's  magnetic 
field  will  be  negligible.  Using  a  potential  gradient  of 
3  mv/km  and  a  crustal  conductivity  of  10  ^  mhos/m,  the 
current  density  is 


—  CT  W  =  (  to'*4)  ^*X|cT  •=  3X(o,0cv~f*//mt5 .2 
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Assuming  that  the  current  forms  an  infinite  sheet  of  20  km 
thickness,  the  magnetic  field  is 

H  =  -L(Xx)£  ^0  —  C5*16  '*)  (s*  =  3aio4!«^.|^4-ua.ws /£►, 

6a  K  £4xio‘*  OuskckC^  AU»S  ) 

This  field  is  completely  negligible.  However,  it  is 
possible  that  the  thermoelectric  currents  are  much  more 
important  inside  the  mantle  due  to  its  high  conductivity, 
but  then  our  assumptions  concerning  the  geothermal  gradient 
and  magnitude  of  the  coupling  coefficients  are  untenable. 

Section  II  -  S.  P.  Currents  Caused  by  Electrochemical 
Coupling 

Electrochemical  action  is  believed  to  be  by  far  the 
most  important  mechanism  of  generation  of  self -potential 
currents  that  are  commonly  associated  with  ore  bodies  or 
other  inhomogeneities  in  rock  properties.  These  currents 
are  assumed  to  arise  from  concentration  differences  of 
electrolytic  solutions  in  contact  with  metallic  conductors 
(such  as  sulfide  ores),  from  the  chemical  differences  of 
the  materials  coming  in  contact  with  solutions  or  from 
differences  in  diffusion  properties  of  materials.  Because 
many  of  the  large  anomalies  that  have  been  observed  are 
associated  with  sulfide  ore  bodies,  nearly  all  the  proposed 
theories  of  the  electrochemical  origin  of  self -potentials 
depend: heavily  on  the  oxidation-reduction  reactions  involv¬ 
ing  sulfides.  Then  theories  are  mainly  concerned  with  the 
mechanism  of  generation  of  self-potentials,  and  few  of  them 
consider  the  effect  of  geometry  on  the  magnitude  of  these 
potentials. 


-  84  - 


In  Chapter  II,  we  studied  the  role  of  electrochemical 
diffusion  in  producing  self -potentials ,  and  in  Chapter  III, 
it  was  shown  that  whenever  inhomogeneities  in  rock  properties 
cut  across  the  surfaces  of  equal  chemical  potential, 
electric  currents  are  generated  through  diffusion.  The 
magnitude  of  the  currents  depends  on  the  differences  in 
coupling  properties  and  the  chemical  potentials  of  the  ions. 

Electrochemical  Coupling  Coefficients  -  Few  measurements 
of  electrochemical  coefficients  can  be  found  in  the  liter¬ 
ature,  but  at  the  present  time  the  Geophysical  Laboratory 
at  M.  I.  T.  is  undertaking  such  measurements  on  rock  samples. 
The  few  results  that  are  presently  available  are  shown  in 
Table  5.1. 


Table  5.1  -  Electrochemical  Coefficients  of  Rock 
Samples  (naCl  Solutions) 


Rock  Type 

Sample 

Electro¬ 

Effective  + 

Investi¬ 

Code 

chemical 

Coefficient 

(mv/decade) 

transference t 

gator 

Sedimentary ■ Rocks 

Y3 

7.8  . 

.57 

Heath 

Colorado  Plateau 

Y18 

4.0 

.53 

II 

Y34 

-2.2 

.48  ■■ 

1! 

Y36 

-5.6 

.45 

II 

Y38 

6.8  • 

.56 

II 

Y4l 

3.6 

.53 

II 

Y46 

5.8 

•  55 

II 

' 

Y52 

8.8 

.57 

II 

Altered  Volcanics 

SA  28 

3.8 

.53 

II 

Arizona 

SA  31 

-3.1 

.47 

II 

SA  34 

8.8 

.57 

II 

SA  35 

9.5 

.58 

!  t 

SA  36 

6.8 

.56 

II 

Quartz-Latite 

BC  4 

6.0 

.55 

It 

Porphyry,  New 

BC  5 

"4.8 

.54 

II 

Mexico 

BC  9 

8.8 

.57 

ti 
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Me tamorphics, 

N  12 

5.8 

.55 

II 

Ontario 

N  16 

5.8 

•  55 

II 

N  37 

6.8 

.56 

It 

Dakota  Sandstone 

S  22A 

-0.7 

.49 

It 

Dense  Sandstone 

G  2 

6.8 

.56 

II 

tuff 

1.2 

.51 

Marshall 

tuff 

26.0 

.72 

tl 

Tremolite 

Sandstone 

43.7 

.87 

II 

Sandstone,  med. 

grained 

-2.4 

.48 

II 

Sandstone,  med. 

grained 

-1.2 

.49 

It 

Sandstone,  fine 

grained 

46.1 

.89 

II 

Sandstone,  dirty 

11.8 

.60 

It 

■The  measurements  were  made  by  placing  the  samples  between 

r* 

two  aqueous  NaCl  solutions  of  strength  .001N  and  .03N,  and 
measuring  the  difference  of  coupling  between  the  rock  sample 
and  water.  The  results  reported  here  are  in  millivolts  per 
decade  of  concentration  (decade  is  the  logarithm  to  the 
base  ten  of  the  ratio  of  concentrations),  and  have  been 
corrected  for  the  coupling  effect  of  water. 

It  is  interesting  to  compare  these  measurements  with 
the  maximum  coupling  coefficient  possible  for  a.  uni -uni valent 
salt  solution.  From  2.46  and  2.48-,  the  coupling  effect  for 
a  two-ion  system  is 


cu  =  huh* iad; 

/  .«*'  F 

Using  2.18  and  the  definition  of  transference  numbers  -fc4 
and  i"  ,  the  diffusion  coefficient,  as  used  in  Table  5.1, 
is  at  25°C  given  by 


5.4 


C  =  25  ZJ  (+tt") 
F 


5.5 


5-^.  \  (  t*-  t")  rrw /decade 


or 


C  ~ 
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The  maximum  value  of  the  coupling  coefficient ,  for  tt  =  1 
and  t  =  0,  is  therefore  59  mv /decade.  Equation  5.5  has 
been  used  to  compute  the  transference  numbers  listed  in 
Table  5-1.  The  largest  reported  values  of  are  close  to 
unity  (forcflays)  corresponding  to  coupling  values  of 
slightly  less  than  59  mv/decade. 

Dependence  of  the  Chemical  Potential  on  Depth  -  The  chemi¬ 
cal  potential  of  an  ion  in  solution  is  essentially  a 
function  of  temperature  pressure,  and  concentration.  Near 
the  surface  of  the  earth,  these  conditions  are  found  to 
be  very  much  dependent  on  depth.  The  temperature  and  the 
pressure  can  be  considered  to  vary  almost  linearly  with 

rs 

depth.  "  The  dependence  of  concentrations  on  depth  is  not 
so  simple,  due  to  the  large  number  of  controlling  factors. 

The  lateral  variations . of  these  factors  are  so  great  that 

?  > 

they  tend  to  mask  the  dependence  of  concentrations  on  depth. 
There  are  very  few  data  on  the  composition  and  chemi- 

A- 

cal  properties  of  pore  fluids  at  depth.  Most  analysis  of 

rocks  give  the  total  ion- content  (solid  +  dissolved)  and 

thus  provide  no  information  on  the  composition  of  the  pore 

fluids.;;  However,  there  are  many  analysis  of  shallow 

ground  ;waters  and  of  oil  field  brines  because  of  their 

obvious  applications.  Table  5.2  lists  a  few  analysis  of 

surface' ground  waters  (Hawkes  and  Webb,  1962)  and  oil 

field  brines  from  Illinois  Basin  (Meents  and  al.,  1952). 

Although  the  data  presented  here  is  very  incomplete,  it 

gives  an  idea  of  the  wide  range  of  variations.  It  is 

obvious  from  these  data  that  the  concentration  of  each  ion 

as  well  as  the  total  salinity  can  vary  by  as  much  as  a  lO-' 

4 

or  10  from  one  condition  to  a  different  one. 

We ; are  here  interested  in  variations  over  depths  of  a 
few  hundred  feet  or  a  few  thousand  feet,  and  for  these 
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depths,  the  most  important  factor  controlling  the  concentra¬ 
tions  is  the  weathering.  It  is  Important  to  realize  that 
it  is  not  necessary  to  know  how  these  changes  are  brought 
about,  in  order  to  apply  the  electrochemical  diffusion 
theory.  All  that  is  required  is  to  have  a  gradient  of 
chemical  potential  that  will  induce  ionic  flows.  As  a 
matter  of  fact,  we  can  assume  that  the  chemical  potential 
varies  very  rapidly  over  a  certain  depth  (the  weathering 
zone)  and  then  varies  more  slowly  or  remains  the  same 
below  that  depth  (depth  environment).  A  study  by  Sato 
(Sato,  1960a;  Sato  and  Mooney,  i960)  of  oxidation  potential 
Eh  and  hydrogen  concentration  pH  of  natural  waters  in 
several  mines  indicates  that  the  Eh  and  pH  vary  over  wide 
ranges  .in  the  weathering  zone  (corresponding  to  changes 
of  concentration  for  Fi***  and  as  high  as  10 and  10^) 
while  they  are  restricted  to  much  narrower  ranges  in  the 
depth  environment.  This  result  would  tend  to  confirm  our 
picture; of  the  depth  dependence  of  the  chemical  potential. 

When  the  system  contains  only  two  kinds  of  ions  of 
equal  valency,  one  positive  and  one  negative,  there  is  only 
one  chqpiical  potential  to  consider  because  the  concentra¬ 
tions  qf  positive  and  negative  ions  can  never  by  very 
different  (  Mn  )•  this  case,  the  coupling 

coefficient  for  tne  medium  is  essentially  constant  and  Is 
given  t»y  5*5*  The  dependence  on  depth  of  the  chemical 
potential  can  be  represented  by 

/  .  -K*  5  6 

5 

where  1  ‘/«  is  related  to  the  depth  of  weathering  (this 
choice  ;pf  z-dependence  is  prompted  by  its  analytic  simplicity) . 
We  have' arbitrarily  assumed  that  at  depth  the  chemical 


Table  5*2  -  Composition  of  Ground  Waters  and  Oil  Field  Brines 
Total  Ions  and  Major  Constituents  (ppm)  (Data  from 
Hawkes  &  Webb,  1962;  Meents  and  al.,  1952) 
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potential  is  zero;  adding  a  constant  value  will  not  affect 
the  induced  currents.  (  represents  the  difference  in 
chemical  potential  and  is  given  approximatively  by 

(Am)  =  RT  JU  Sill3  5.7 

/  Cto) 

where  and  C Co)  art  the  concentrations  at  depth  and  at 

the  surface  respectively.  We  have  here  neglected  the 
variations  of  jU  due  to  temperature  and  pressure  changes. 
Subject  to  these  restrictions,  the  equations  derived  in 
Section  III  of  Chapter  III  can  be  used  to  compute  the 
electrochemical  diffusion  potentials  of  solutions  contain¬ 
ing  only  two  kinds  of  ions  (or  when  two  ions  predominate). 

When  several  ions  are.  present  in  appreciable  quantities, 
the  concentrations  are  not  necessarily  equal,  and  we  can 
not  speak  anymore  of  the  variations  of  a  single  chemical 
potential.  Starting  from  equation  2.53>  giving  the  current 
flow  in  a  multi-ion  system,  the  potential  is  defined 
now  as' 


{  =  /  +  5.8 

The  electric  potential  p/  and  the  chemical  potentials  jj. • 
must  be  continuous  across  any  boundary,  therefore,  using 
2.18  and  2.20,  we  get  for  the  boundary  condition  on  -f 


U-C>£a- 

i 


where 
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Uj t  C'^and  yU,'  are  the  mobility,  the  concentration,  and 
the  chemical  potential  of  each  ion,  respectively;  sub¬ 
scripts  and  superscripts  1  and  2  refers  to  the  two  media 
in  contact.  This  equation  replaces  the  boundary  equation 
2.58. 

The  above  expression  cannot  be  separated  into  the 
product  of  a  coupling  difference  and  a  driving  potential 
as  we  did  in  the  case  of  a  dissociated  salt.  In 
principle  we  can  compute  its  value  from  a  knowledge  of 
the  mobilities  and  the  concentrations  of  all  the  ions 
present.  Since  the  mobilities  and  the  concentrations  depend 
on  depth,  we  can  still  assume  that  the  depth-dependence  of 
the  expression  on  the  right  hand  side  of  5.9  is  given  by 
5.6. 

Electrochemical  Diffusion  Potential  across  a  Dike  -  In  the 
absence  of  metallic  conductors,  the  electric  potential  is 
due  to  diffusion  of  ions,  very  much  similar  to  a  junction 
potential.  But  when  metallic  conductors  such  as  sulfide 
ore  bodies  are  present,  some  complications  are  Introduced 
which  will  be  discussed  in  the  next  section. 

The  general  features  of  the  electrochemical  diffusion 
potential  across  a  dike  were  computed  in  Chapter  III,  and 
are  shown  in  Figs.  3.12  and  3.13.  The  parameter  /K 
represents  roughly  the  depth  of  weathering  or  the  depth 
associated  with  large  changes  in  concentrations;  therefore, 
the  product  (  K  d  )  represents  the  ratio  of  the  dike 
thickness  d  to  the  weathering  depth.  In  Fig.  3*12  the 
geometric  factor  <f>(x)  has  been  plotted  as 

a  function  of  (jj  where  x  is  the  distance  along  the 
surface.  These  plots  show  a  great  similarity  in  amplitude 
and  the  rate  of  decrease  with  distance.  Clearly,  the 
maximum  value  of  the  geometric  factor  is  only  slightly 
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dependent  upon  the  parameter  (  Kd  )  of  the  dike  (i.e.  the 
thickness  of  the  dike).  In  going  from  a  thick  dike,(  Kd  )= 
2.0,  to  a  thin  dike,  (W)  =  0.5,  the  maximum  value  of  the 
potential  changes  only  by  a  factor  of  two. 

The  dependence  of  the  potentials  on  conductivity 
contrast  fit|  is  much  more  important,  especially  for  thin 
dikes.  For  a  dike  of  thickness  (  Kci  )  =  0.13,  the  maximum 
potential  changes  by  a  factor  of  four  in  going  from  a 
conductive  dike,  fin  =  0.5*  to  a  more  resistive  dike,  R^.-o.r. 
The  changes  are  even  more  dramatic  for  larger  conductivity 
contrasts.  For  a  highly  resistive  dike  (  =-l-o  ), 

the  electrochemical  coupling  effect  is  almost  negligeable 
as  it  flan  be  seen  from  equations  3.37a  and  3* 37b,  while 
for  a  very  conductive  dike:  (  ^  +  i.o  ),  the  maximum 

value  of  the  geometric  factor  is  almost  equal  to  unity 
and  is  not  affected  appreciably  by  the  thickness  of  the 
dike.  "This  is  one  of  the  reasons  for  the  larger  potentials 
associated  with  conductive  ore  veins. 

Figure  3.12  shows  that  the  half -power  width  of  the 
anomaly,  for  thin  dikes  computed  cases,  is  almost  independent 
6f  the  "parameter  (  Kel  ),  implying  that  for  these  cases,  the 
half -power  width  depends  on  the  thickness  of  the  dike. 

For  very  thick  dikes  the  vertical  contact  case  (Fig.  3.11) 
is  applicable . 

The  magnitude  of  the  electrochemical  potential  can  be 
computed  from  a  knowledge  of  the  mobilities  and  the 
concentrations  of  the  ions  present.  Unfortunately,  the 
pore  fluids  of  rocks  represent  multi-ion  systems  whose 
coupling  properties  have  not  been  measured.  The  coupling 
coefficients  listed  in  Table  5*1  have  been  measured  using 
a  NaCI ' solution. 

As,  a  first  approximation,  we’  can  assume  that  only  two 
kinds  of  ions  (e.g.  Na+  and  Cl”)  are  present.  This  is 
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very  often  a  good  approximation  as  can  be  seen  from  the 
data  in  Table  5.2.  Therefore,  we  can  use  equations  5.6 
and  5-7  which  are  valid  for  a  two-ion  system.  The 
concentration  C  in  5-7  would  refer,  for  instance,  to 
the  concentration  of  sodium  or  chlorine  or  the  total 
concentration.  From  Table  5-2,  it  is  seen  that  the  total 

P 

ion  concentration  can  vary  easily  by  a  factor  of  10  over 
a  depth  of  two  or  three  hundred  feet,  giving  rise  to  a 
driving  potential  b.jA  -  RTAwo2  .  Corresponding  to  this 
choice  of  concentration  difference,  the  depth  parameter 
!/k  of  equation  5-6  would  be  about  100  ft.  or  more. 

From  Table  5*1j  typical  values  of  the  coupling  differ¬ 
ence  (C,-C,,  )  are  seen  to  be  about  5  to  15  mv/decade. 
Assuming  a  dike  thickness  of  100  ft.  (  i.o  ) }  a 
conductivity  contrast  of  three  (  °--5"  ),  and  a 

coupling  difference  of  15  mv/decade,  we  get  for  the  maxi¬ 
mum  diffusion  potential  (see  Fig.  3.12) 

A</>  =.  (iff-  m'j/dtcaJj  )(Ug  |ol)(o.t)  s  IS  mV  5-10 

This  can  be  considered  as  a  typical  maximum  value  of 
the  potentials  induced  by  electrochemical  coupling.  It  is 
rather  'difficult  to  produce  much  larger  potentials  by 
diffusion. 

It K  can  be  seen  from  5 >9  that  the  driving  potential  is 
a  weighted  average  of  all  s  and  the  weighting  factors 
are  proportional  to  the  concentrations  of  the  ions.  If  the 
concentration  of  an  ion  goes  to  zero,  it  is  true  that  its 
chemical  potential  becomes  infinitely  large,  but  the  product 
in  equation  5.9  goes  to  zero.  This  point  is  very 
important  in  understanding  why  huge  potential  differences 
are  never  produced  by  diffusion  coupling. 
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Role  of  Metallic  Conductors  -  Most  of  the  observed  self- 
potentials  are  less  than  100  mv  in  magnitude.  This  certainly 
is  true  of  the  background  noise  in  S.  P.  surveys.  However, 
larger  potentials  of  several  hundred  millivolts  are  often 
observed  over  sulfide  ores  and  graphite  zones  that  form 
massive  bodies  or  continuous  veins.  These  large  self- 
potentials,  that  can  attain  500  mv  in  magnitude,  can  not 
be  caused  by  electrochemical  diffusion  alone.  In  the 
preceeding  section,  we  found  that  the  largest  anomalies 
induced  by  simple  diffusion  can  not  be  more  than  tens  of 
millivolts  unless  rather  unreasonably  high  value  of 
are  assumed.  When  metallic  conductors  are  present,  such 
an  increase  in  the  driving; potential  is  produced. 

We.  shall  see  farther,  that  in  the  presence  of  electronic 
conductors,  the  driving  potential  is  essentially  determined 
by  the  concentrations  of  the  minor  elements  capable  of 
undergoing  oxidation-reduction  reactions,  such  as  iron  and 
manganese.  Among  these  minor  constituents,  iron  occurs  in 
rocks  in  fair  quantities,  however,  its  solubility  is  very 
low.  Prom  a  knowledge  of  all  the  reactions  and  the  free 
energies  Involving  iron  we  can  compute  its  ionic  concen¬ 
tration.  These  well-known  methods  of  chemical  thermodynamics 
need  not  be  repeated  here,  and.,  the  Interested  reader  is 
referred  to  R.  H,  Garrels  (I960)  and  M.  J.  N.  Pourbaix  (19^9) 
for  detailed  treatments.  We  shall  confine  ourselves  to  a 
brief  discussion  of  the  results. 

It  is  found  that  any  chemical  reaction  in  aqueous 
solutions  can  be  represented  in  terms  of  the  reacting  species, 
the  pH,  and  the  Eh  of  the  solution .*Ji* pH,  defined  jas  -log  (H+), 
determines  the  influence  of  the 
oxidation  potential,  determines 
For  instance,  the  activities  of 
equilibrium  with  each  other  are 
Mooney;,  I960,  p.  232) 

R*-"  ^=4  P.”*+  4- 

01 -"7I 


hydrogen  ion,  and* Eh,  the 
the  influence  of  electrons . 
ferrous  and  ferric  ions  in 
determined  by  (Sato  and 

r. to- 1 
5".  io.  Z 
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Therefore,  a  change  of  0.12  volts  In  the  Eh  value  changes 
the  ratio  of  ferric  to  ferrous  ions  by  two  orders  of 
magnitude. 

If  the  equation  5.10.1  represents  the  only  important 
reaction,  the  concentrations  of  ferric  and  ferrous  ions 
are  independent  of  the  pH.  But,  in  general,  there  are 
other  competing  reactions  such  as 

+  5.10.3 

3PH  5-10.4 

Equation  5*10.4  implies  that  at  constant  Eh,  a  change  of 
unity  in  the  value  of  the  pH  corresponds  to  a  change  of 

o  1 

10“*  in  the  concentration  of  ions.  Actually  because 

,  of  the  many  competing  reactions  the  concentrations  of 
ferrous  and  ferric  ions  are  determined  by  the  combination 
of  both  pH.  It  is  found  that  in  very  oxidizing  environ¬ 
ment  values  of  Eh  =  0.7  volts  and  pH  =>  2  are  possible; 
while,  in  very  reducing  conditions  values  of  Eh  =  -0.3 
volts  and  pH  =  8  are  possible  .  Using  these  values  of 
Eh  and  pH  in  conjunction  with  the  data  presented  by 
R.  M.  Garrels  (Garrels,  I960,  Pigs.  6.5a  and  6-.5d),  we 
see  that  the  concentrations  of  Fe*5  and  Ft*2,  vary  respec¬ 
tively  from  10  ^  and  10  ^  at  the  surface  to  10  2  and  10  ^ 
at  depth.  When  electronic  conductors  such  as  sulfide  ores 
or  graphite  are  present  it  will  be  shown  that  the  magnitude 
of  the  S.  P.  currents  are  controlled  by  the  chemical  poten¬ 
tials  of  the  reacting  species  (i.e.  fv* 3  and  )  and  with 
such  large  variations  in  the  concentration,  large  S.  P. 
anomalies  are  possible. 

Many  sulfide  ores  such  as  pyrite  and  chaleophrite  are 
good  electronic  conductors,  and  when  present  in  a  massif  form 
or  as  continuous  veins,  they  provide  a  conducting  path 
for  elections  between  the  surface  where  oxidizing  conditions 
prevail  and  the  zone  below .the  water  table  where  reducing 
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conditions  are  prevailing.  The  presence  of  this  conduct¬ 
ing  path  facilitates  the  oxidation- reduction  reactions. 

Let  us  specifically  consider  the  ferric-ferrous  ion  pair. 
Below  the  water  table,  the  reaction  taking  place  is 
Ft** — s>  +  e~  which  gives  off  an  electron  that 

travels  up  through  the  ore  body.  In  the  weathering  zone 
above  the  water  table  the  opposite  reaction  takes  place: 

Fi+  3-i~Q  _ *  .  As  a  net  result  of  then  two 

reactions,  a  ion  disappears  at  depth  and  reappears 

near  the  surface,  while  a  Ft'*''3  ion  disappears  from  near 
the  surface  and.  reappears  at  depth.  This  gives  the  fei"h 
and  fit*3  ions  apparent  mobilities  that  are  much  higher 
than  the  mobilities  of  the  diffusion  ions .  We  shall 
designate  these  apparent  mobilities  as  11*  .  Of  course, 

inside  the  ore  body  itself  the  other  ions  have  zero  mobility. 

Iff.  other  reacting  species  are  present  (e.g.  Manganese) , 
they  behave  in  very  much  the  same  way  as  the  pair. 

Sato  and  Mooney  (i960)  have  studied  the  chemical  reactions 
taking iplace  around  sulfide  ore  bodies,  and  their  study 
confirm^  the  importance  of  the  ferrous-ferric  electron 
exchange  reaction. 

According  to  equation  5.9,  the  driving  force  is  a 
weighted  average  of  the  chemical  potentials  yUy  of  all  the 
ions  present,  the  weighting  factors  being  dependent  on  the 
mobilities.  Inside  the  conducting  ore  body,  only  the 
reacting  species  have  large  apparent  mobilities  and  all 
other  ions  have  zero  mobility  (it  is  obvious  that  the  ore 
body  does  not  influence,  for  instance,  the  diffusion  of C^f). 
Therefore,  if  the  ferrous-ferric  pair  represents  the  only 
important  reacting  species,  the  driving  force  ~YX  in  the 
ore  zone  (medium^)  is 


5.11a 
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where 


P  +  Mst*^**^** j  5 -lib 

r$taM  t/tmwi/njt,  t/m/a&  ao^  The  conductivity  seems  to 

be  dependent  on  the  concentrations  of  ferrous  and  ferric 
ion.  This  is  an  apparent  inconsistency  that  can  be  easily 
removed.  For  one  thing,  we  know  that  inside  the  ore  body, 
the  conductivity  is  electronic  and  should  not  depend  on  the 
ferrous-ferric  pair.  Furthermore,  while  the  Hr3  and  Ftft- 
ions  have  essentially  infinite  mobilities  through  the  ore 
body  (we  are  neglecting  the  reaction  impedance),  their 
concentrations  are  zero,  which  implies  that  the  products 
and  are  constant  and  independent  of  the 

mobilities  and  the  concentrations.  Finally  the  ferrous 
and  ferric  ions  are  coupled  through  the  reaction  ^ /5+Ve_ 
for  every  ferric  ion  transferred,  a  ferrous  ion  is  trans¬ 
ferred  in  the  opposite  direction — which  means  that  the 
mobilities  and  concentrations  inside  the  ore  are  related 
by 


-  -  uE*  cu 


F'+'^Fe*' 


5.12 


Substituting  from  5.12  and  5.11b  into  5.11a,  we  get 


f  1  iv] 


5.13a 


(in  millivolts)  5.13b 


This  above  expression  represents  the  driving  potential 
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inside  the  ore  body  (medium  2)  when  the  ferrous-ferric 
pair  represents  the  only  important  reacting  species. 


Outside  the  ore  body,  the  driving  potential  is  still 
given  by  5.9 


y  ,  y  Ut  Ci  Mi 

71  t  ^ 


5.14 


where  all  the  ions  have  finite  mobilities.  If  the  concen¬ 
tration  of  a  particular  ion  is  very  small,  its  contribution 
to  the  driving  potential  will  be  equally  small  because 

C:  Li i' - *>o  .  The  .  magnitude  of  the  driving  potential 

/  •:  c,4o 

is  determined  mainly  by  the  ions  that  are  most  abundant  . 
(e.g.  ns,  cr,cs<-  etc.),  and  the  concentrations  of  these 
ions  vary  only  slightly  with  depth.  From  the  data  in 
Table  ^.2,  it  is  obvious  that  at  best  these  concentrations 
vary  by  two  or  three  orders  of  magnitudes  which  can  only 
produce  small  anomalies.  Therefore,  the  variations  of 
with  depth  are  essentially  negligible,  and  we  can  write 
for  the  total  driving  potential 


K 


5~%  I 


5.15 


This  last  equation  is  similar  to  the  Nerst's  equation 
for  reversible  cell  reactions,  however,  the  two  situations 
are  not  quite  equivalent.  In  the  reversible  cell  the  con¬ 
centrations  are  constant  throughout  the  cell  and  the  e.m.f. 
is ' determined  by  the  equilibrium  between  ferrous  and  ferric 
ions.  In  our  case,  equation  5.15  represents  a  variable 
boundary  condition  as  both  fzt3  and  fy^are  functions  of 
depth. 
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It  is  by  now  clear  why  such  minor  constituents  as 
iron  that  are  normally  unimportant  in  determining  the 
electrochemical  properties  of  the  soil,  suddenly  become 
so  important  in  the  presence  of  a  metallic  ore  body. 

Another  feature  of  this  explanation  is  that  it  does  not 
require  the  ore  body  to  undergo  oxidation  or  reduction. 

The  ore  body  provides  only  a  conducting  path  for  the 
electrons.  This  last  conclusion  is  consistent  with  the 
arguments  advanced  by  Sato  and  Mooney., 

The  equation  5.15  is  in  the  form  of  the  product  of 
a  coupling  difference  (  c.(-Ct  =  59.1  mv/decade)  and  a 
driving  potential  (  Ayu  =  log  )•  Therefore,  we 

can  sti'-ll  use  the  '  equations  3.28,  3.27a,  3.27b,  and  the 
Pig.  3.-12. 

Prom  a  knowledge  of  the  Eh  and  pH  values  around  the 

ore  body,  we  can  obtain  the  concentrations  of  ferrous  and 

ferric rions,  and  using  5.15  and  Fig.  3.12,  we  can  compute 

the  S .  ;P .  anomalies  to  be  expected.  Values  of  Eh  =  0.6  and 

pH  =  4  ;are  typical  for  the  near  surface  conditions  (Sato 

and  Mooney,  I960;  Fig.  4),  while  values  of  Eh  =  -0.2 

and  pH  -  8  are  typical  of  the  conditions  at  depth.  Using 

these  values  of  Eh  and  pH,  it  is  found  (Garrels,  I960; 

Figs.  6.5a,  6.5d,  and  6.6)  that  in  going  from  the  oxidizing 

conditions  near  the  surface  -to  the  reducing  conditions  at 

depth,  the  ferric  ion  concentration  decreases  approximatively 
12 

by  10  while  the  ferrous  ion  concentration  increases  by 
2 

about  10  .  Therefore,  the  driving  potential  difference  is 
K  1  lO_,#  1  ,  -‘I 

V  =  l’5^=  *”5'“ 

Assuming  a  weathering  depth  of  500  ft.  (  V4c  =  and  a 

conductive  dike  <at  250  ft.  thick  (  x4-0.iT  ), 
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the  geometric  factor  is  seen  to  be  about  0.6.  Using  these 
values,  the  maximum  S.  P.  anomaly  observed  over  the  ore 
body  is 

A</=  (ST  00*3  io'z)(o.O  =^25' 

This  certainly  is  a  very  large  self -potential  and 
compares  favorably  with  some  of  the  largest  observed  S.  P. 
anomalies.  Most  S.  P.  anomalies  are  below  500  mv,  but 
occasionally  larger  S.  P.  anomalies  have  been  observed. 

Sato  mentions  some  unusually  large  values  reported  by 
Shibato  (1.0  volt,  graphite)  and  Fujita  (1.3  volts). 

Assuming  extremely  oxidizing  reducing  conditions  near 

the  surface  (Eh  =  0.7  volts,  pH  =  2.5)  and  extremely 

reducing  conditions  at  depth  (Eh  =  -0.3  volts,  pH  =  9), 

the  ferric  ion  concentration  decreases  by  lO1^  and  the 

•5 

ferrous  ion  concentration  decreases  by  10J  in  going  from 
the  oxidizing  zone  to  the  depth  environment.  If  the  ore 
body  is  extremely  conductive  (  **  1.0),  the  geometric 

factor  ;is  near  unity  (see  equation  3 -37a  which  for  (7,,-  1.0 
and  X  4,0,  gives  f,  (x)  =  (c,~  CiXfy))-  Therefore,  the  maximum 
anomaly  is 

j£X'-<0  -  <,so 


Thp.s  is  just  about  the  largest  anomaly  we  can  get 
through  the  oxidation-reduction  reactions  of  the  ferrous- 
ferric  ;lon  pair.  If  the  larger  observed  S.  P.  anomalies 
are  controlled  by  oxidation-reduction  reactions,  other  ions 
must  be  involved.  Mh  N*tZ  pair  is  one  possibility,  however, 
it  is  not  clear  under  what  conditions,  the  manganese  ions 
become;,  important.  This  point  needs  to  be  investigated 
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further.  It  is  also  possible  that  some  of  the  larger 
anomalies  are  combinations  of  electrochemical  S.  P. 
potentials  and  topographic  effects  due  to  electrofiltra¬ 
tion. 

Section  III  -  Effect  of  Random  Inhomogeneities  on  S .  P. 

Currents 

Although  the  diffusion  coupling  properties  of  rocks 
do  not  seem  capable  of  producing  the  large  S.  P.  anomalies 
observed  over  certain  ore  zones.  They  can  still  give  rise 
to  some  S.  P.  variations  which  must  be  considered  as  noise 
when  S.  P.  surveys  are  being  conducted.  In  this  section, 
we  shall  evaluate  the  effect  of  a  random  distribution  of 
inhomogeneities  in  producing  such  S.  P.  noise,  when  the 
variations  in  coupling  properties  are  function  of  the 
x-coordinate  only. 

__  The  general  expression  for  the  rms  value  of  the  noise 
fj"  was  derived  in  Section  IV  of  Chapter  III.  From  3.49, 
we  have 


ft  =  "IPX  i 

5.17 

where  ^.p(W'rl)  is  the  autocorrelation  of  the  potential  due 
to  a  vertical  contact;  ^0T0  is  the  autocorrelation  of 
the  gradient  c<  (x )  of  the  coupling  parameter  C 
and  I/k-  is  essentially  the  depth  of  weathering. 

The  particularly  simple  model,  that  we  shall  use  to 
compute  OL(l'cl)  is  that  of  an  earth  composed  of  an  infinite 

*  i 

number  of  dikes,  each  having  a  thickness  t:  and  a  coupling 

A 
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coefficient  C-  that  we  only  know  by  their  statistical 

A tr 

properties . 

As  we  cross  each  dike  boundary,  the  coupling  coeffi¬ 
cient  changes  by  a  finite  amount  AC;  which  can  be  positive 
or  negative.  Thus,  the  gradient  *(x )  of  coupling  is 
represented  by  a  series  of  impulses.  The  variation  of  o^(x) 
across  a  dike  is  shown  in  Pig.  5.1.  The  average  value  of 
the  coupling  difference  AC;  is  obviously  zero,  but  its 
rms  value  is  not.  We  theft  have 

=  o 

»  (*z)Z  5,18 

where  (kd)  is  the  rms  value  of  the  fluctuations  in  coupling 
properties . 


Figure  5.1  -  Typical  variations  of  coupling  gradient  c/(x) 

across  a  dike. 


The  distribution  of  the  random  pulses,  along  the  surface 
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depends  on  the  random  distribution  of  the  dikes.  From  a 
well-known  theorem  in  communication  theory,  the  random  • 
function  «A-(x)  ce  ie  represented  as  a  convolution  of  a 
transient  function  sl(x)  and  white  noise  function 

OQ 

(*)  =  j  I  (x- u}  du  5.19 


We  shall  assume  that  in  our  model,  A(.x)  can  be  represented 
by  anology  with  Fig.  5-1,  by  an  impulse  function 

/lLx)  =  (&M)  [  &(x)  5.20 

where  ~t  is  related  to  the  average  dike  thickness  and  (AM) 
is  related  to  the  rms  fluctuations  in  coupling  properties 
(AMC  =  o)  . 

With  these  assumptions,  the  autocorrelation  is 

given  by 


Cx)  <xx  r) 


Jj  ^0*9  •’tCn*)  y^ctx- u'J  t-M'9  ciu’^w  w 


00 

-  jj  >L(a’)  AU")  [  L,  S(*-u-r)J  duct*  " 

-art  0 


1  j  A(tt’)  fc(utc)d 


°r  'L(Iti)  =  5.21 

~h  indicates  the  number  of  dikes  per  unit  length,  and  (t) 

t>  y-f- 

is  the  autocorrelation  of  tt-Cx)  .  From  5*20,  we  get  for  <§r(x) 


5rr(r')  =  C^M")2  [iS(r)- 


5.22 


and  substituting  from  5.21  and  5.22  into  5.17>  we  get  for 
the  S.  P.  noise 

[few-H  5.a3 

The  dimensionless  parameter  (Kt)  in  the  above  equation  is 
a  measure  of  the  ratio  of  the  dike  thickness  to  the  depth 
of  weathering.  We  finally  have  to  relate  (AM)*  to  the  rms 
fluctuations  of  the  coupling. 

If  equation  5.20,  giving  *.(.*) ,  represents  the  typical 
change  of  the  coupling  gradient  M  ,  the  typical  change 
in  the  .coupling  itself .  Would  be 

i 

A  CA-  =  (am)V(O  5.24 

V(x)a 

whereAisApulse  of  finite  width  "t  and  height  of  unity. 
Therefore,  the  variations  of  A along  the  surface  can 
be  represented  in  our  model  by 

=  j  (AM)  5.25 


and  its  autocorrelation  is  given  by 

O0 

A C;6o  Ac/x-rr)  =  jj  (too'UhOVCv)  skO<-u')A;lxtT-^)  Jl 


u  'du 
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The  rms  value  of  fluctuations  is  obtained  by  putting  =  © 
(Acf=  A  C]  =  5.26 

Substituting  from  5.26  into  5.23,  we  get  the  final  expression 


For  very  thick  dikes  the  numerator  in  the  brackets 

remains  finite  while  the  denominator  becomes  infinitely 
large,  and  thus,  the  rms  valye  of  the  noise  approaches 
zero.  At  the  other  end,  when  the  dike  becomes  infinitely 
thin  (  kt-»o  )  the  noise  level  ft  approaches  a  limiting 
value  given  by 

5-28 

where  is  the  derivative  of  the  autocorrelation 

function  shown  in  Fig.  3.14.  In  this  case,  the  rms 
fluctuations  (Ac )*  of  the  coupling  properties  becomes 
small  in  amplitude,  since  the  adjacent  dikes  cannot  be 
very  different  in  their  properties.  The  noise  will 
consists  essentially  of  high  frequency  components.  For 
large  values  of  the  thickness  "b  ,  the  same  high  frequency 
components  may  be  present,  but  there  are  also  some  broad 
variations  in  coupling  properties  which  will  give  rise 
to  low  frequency  components  in  the  noise . 

USing  a  coupling  difference  of  10  mv/decade,  a  con- 
centration  ratio  of  10  ,  and  a  thin  dike  ,  we  get 

f  ■  &  Ugio*-/  r  2.g-o.q  1 
x  y|“»-  l  bo  -* 
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and  a  rms  value 


i 


9  mv 


Although  this  value  of  19  mv  should  be  considered  as  the 
upper  bound  of  typical  diffusion  noise  to  be  expected,  it 
is  large  enough  to  account  for  much  of  the  observed  noise. 

Section  IV  -  Other  Causes  of  S .  P. 

Electrofiltration  and  bioelectrical  currents  are  also 

believed  to  be  important  in  producing  self-potentials . 

For  instance,  the  topographic  effect  observed  in  S .  P. 

surveys'  is  believed  to  be  due  to  the  streaming  potential  of 

water  seeping  into  the  ground.  We  already  have  seen  that 

rocks  have  large  electrokinetlc  coupling  coefficients  of 

several  millivolts  per  atmosphere.  Therefore,  flow  of 

water  under  a  small  pressure  head  is  sufficient  to  induce 

measurable  potential  differences.  Cagnlard  (1956)  reports 

the  results  of  a  S .  P.  survey  in  Yugoslavia  over  a  very 

porous  quartz  hill,  where  an  unusually  large  observed 

anomaly- of  -300  mv  was  attributed  to  electrofiltration. 

So.me  authors  have  attributed  the  topographic  effect  to 

uphill  currents  caused  by  the  atmospheric  electric  gradient. 

Although  there  is  a  large  vertical  gradient  of  about  100 

volts/n> ,  these  currents  are  negligible  due  to  the  extremely 

-14 

low  conductivity  of  the  air.  Using  a  typical  value. of  10 

mhos/m  for  the  conductivity  of  the  air,  the  Induced  current 
“12  2 

density'  is  10  amps/m  ,  and  the  ohmic  potential  drop 

_4 

inside  [the  ground  with  an  average  conductivity  of  10 
mhos/m  is  therefore  10~^  volts/m  or  0.1  mv  per  100  meters, 
a  potential  drop  too  small  to  affect  S.  P.  results. 

Streaming  potentials  induced  by  flow  of  mud  under 
pressure  are  also  responsible  for  the  self-potentials 
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observed  in  well  loggings. 

Currents  of  bioelectric  origin  are  also  believed  to  be 
very  important  in  producing  large  &.  E.  noise  in  wooded 
areas.  Scott  (1962)  reports  potential  differences  of  20  mv 
or  more  across  plant  roots.  A  series  of  field  studies  of 
S.  P.  noise  by  Madden,  Cantwell,  Sill,  and  Maher  (Madden 
and  Cantwell,  1962)  showed  that  the  micro-structure  of 
the  self -potentials  (noise  level)  was  higher  in  the  woods 
than  in  open  areas.  Several  traverses  run  from  an  open 
field  into  the  adjacent  woods  gave  large  negative  anomalies 
that  were  greatest  right  at  the  edge  of  the  woods  (50 — 100 
mv) .  Bioelectric  currents  are  believed  to  travel  down  the 
main  root  and  return  through  the  earth,  &s  one  would  expect 
for  any -conductor  extending  from  the  surface  to  any  depth. 

Section  V  -  Conclusions 

The  results  of  this  chapter  are  in  general  agreement 
with  the,  observed  facts  as  well  as  the  other  proposed 
mechanisms  of  S .  P.  currents  (Helland,  1946;  Sato  and  Mooney, 
i960;  Wyllie,  1955)-  Many  of  the  results  we  have  obtained 
are  similar  to  the  results  obtained  by  Sato  and  Mooney 
using  th,e  electrochemical  approach.  However,  our  approach 
through 1 >the  use  of  irreversible  thermodynamics  is  superior 
in  certain  respects. 

Fi^st  of  all,  using  Onsager's  relations  as  our  start¬ 
ing  point,  we  were  able  to  include  the  effect  of  geometry 
in  our  equations.  It  is  generally  assumed  that  the  maximum 
anomaly •to  be  expected  in  a  self-potential  survey  would  be 
one  half  of  the  total  potential  difference.  The  computed 
geometric  factors  for  a  dike  (Fig.  3.12)  show  that  this 
rule  of : thumb  is  not  necessarily  true  and  the  actual  magni¬ 
tudes  depend  very  much  on  the  ratio  of  the  dike  thickness 
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to  the  weathering  depth  and  on  the  conductivity  contrast 
For  very  thin  dikes,  the  observed  anomaly  is  only  two  to 
three  tenth  of  the  total  potential  difference,  while  for 
very  thick  and  highly  conductive  dikes  the  observed  anomaly 
is  nearly  equal  to  the  total  potential  difference  available. 

Furthermore,  our  approach  shows  that  S.  P.  anomalies 
can  be  also  caused  by  pure  electrochemical  diffusion.  These 
diffusion  potentials  are  small  compared  to  the  anomalies 
observed  over  certain  sulfide  or  graphite  bodies,  however, 
they  are  of  sufficient  magnitude  to  account  for  much  of 
the  background  S.  P.  noise. 

Finally,  our  treatment  includes  the  ore  bodies  with 
electronic  conductivity  as  an  /special  case.  By  using  the 
irreversible  thermodynamics  approach,  we  have  introduced 
the  concentrations  of  the  reacting  species  (i.e.  the  Eh 
and  the  pH  of  the  solution)  rigorously.  It  is  true  that 
the  Eh-.of  a  solution  containing,  for  instance,  ferric  and 
ferrous  ions  is  determined  by  the  concentrations  of  these 
two  ions,  but  in  the  approach  through  classical  chemistry 
it  is  pot  quite  clear  why  the  Eh  is  related  to  the  observed 
anomalies.  Our  approach,  not  only  shows  clearly  the 
relationship  between  the  activities  of  the  ions  in  solution 
and  the  observed  potentials,  but  it  clearly  shows  the 
limitations.  Briefly  these  limitations  are:  (a)  the  ore 
body  must  provide  a  continuous  path  for  electrons;  (b)  the 
coupling  coefficient  of  59-1  mv/decade  is  only  valid  when 
large  qoncentration  variations  of  one  reacting  species  are 
involved....  When  several  reacting  species  are  present, 
equation  5.15  is  no  longer  valid,  and  new  equation  includ¬ 
ing  all  the  reacting  species  has  to  be  derived. 


CHAPTER  VI 


SUGGESTIONS  FOR  FUTURE  WORK 

In  this  thesis,  the  principles  of  irreversible  thermo¬ 
dynamics  were 'applied  successfully  to  a  study  of  the  broad 
features  of  S .  P.  anomalies  associated  with  underground 
nuclear  detonations  and  sulfide  ore  bodies,  however,  several 
difficulties  were  encountered  that  need  to  be  investigated 
further . 

Interpretation  of  the  Underground  Detonation  Data  -  The 
available  field  data  contained  larger  amount  of  noise  that 
made  the  depth  interpretation  difficult.  As  a  matter  of 
fact,  fair  fits  were  obtained  using  the  curves  computed 
for  several  (j)  ratios,  and  it  looked  as  if  the  data  would 
fit  the  theoretical  curves  for  other  thickness  to  depth 
ratios.  Furthermore,  the  depth  estimates  were  systemati¬ 
cally  too  high.  This  problem  needs  to  be  further  studied 
along  two  different  lines.  First  of  all,  better  field 
recording  techniques  are  needed  to  minimize  the  effect  of 
S.  P.  noise  and  tellurics.  Secondly,  theoretical  curves 
for  many  more  thickness  to  depth  ratios  (  £  )  and  conduc¬ 
tivity  contrasts  need  to  be  computed  to  evaluate 

unambiguously  the  usefulness  of  the  theoretical  model  in 
predicting  the  depth  of  the  shot. 

Chemical  Properties  of  Pore  Fluids  -  Very  little  is  known 
of  the  actual  chemical  properties  of  pore  fluids.  The 
estimates  of  concentrations  used  in  computing  the  S.  P. 
anomalies  over  the  ore  bodies  were  obtained  from  simple 
theoretical  consideration.  The  actual  concentrations  of 
the  ions  in  the  pore  fluids  may  be  very  different  from 
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the  theoretical  values  because  many  reacting  species  may  be 
involved.  An  investigation  of  the  chemical  properties  of 
pore  fluids  is  required  to  obtain  actual  estimates  of  the 
concentrations  of  the  reacting  species.  Occasionally 
anomalies  of  one  volt  or  better  have  been  observed,  but 
these  can  not  be  explained  on  the  basis  of  the  available 
data  on  the  solubility  of  ferric  and  ferrous  ions.  However, 
they  can  be  explained  if  the  pair  is  the  dominat¬ 

ing  reacting  species.  The  solubility  of  manganese  in  pore 
fluids  in  the  presence  of  iron  and  other  ions  needs  to  be 
investigated. 

Multi -Component  Oxidation-Reduction  Reactions  -  In  our  treat¬ 
ment  of ' the  self-potentials  over  metallic  ore  bodies,  we 
assumed  that  ferrous -ferric  ion  pair  is  the  only  important 
reacting  species,  and  the  equation  5-15  that  was  derived 
to  compute  the  potentials  is  only  valid  under  these  condi¬ 
tions.  'The  form  of  this  equation  for  a  system  containing 
several  reacting  species  needs  to  be  Investigated. 

Extension  to  Disseminated  Ore  Bodies  -  In  our  treatment  of 
the  self -potentials  over  metallic  ore  bodies,  we  assumed 
that  the  ore,  was  present  in  massif form  or  as  a  continuous 
vein.  However,  in  most  cases  of  interest  the  ore  is  pre¬ 
sent  in  disseminated  form.  The  reacting  species  may  still 
have  large  apparent  mobilities  and  dominate  over  the  other 
ions.  This  problem  needs  to  be  investigated  because  of 
its  practical  applications. 
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Solution  of  the  Boundary  Value  Cases 


I .  One  Layer  over  a  Half -Space. 
Boundary  conditions  are 
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II.  Two  layer  over  a  Half -Space, 
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Solving  these  equations,  we  get 
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III,  Vertical  Contact. 
We  have 


L  .  r 

'  v  ~  o 


«x>-. 


X<t> 


where  ft  and  ft 


satisfy  the  boundary  conditions 


-  113  - 


'bf. 

.  *9. 

-  — -*  -=>■  o 

3-  ss.  i 

Tsi 

~6b 

*£l 

a  <J-  5i«- 

tx 

tx  | 

1  <x}x^  t> 

Applying  the  first  boundary  condition  we  get 
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and  a  similar  equation  for  Tx  .  Therefore  we  get 
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Applying  the  other  boundary  conditions  we  get 
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We  also  have  the  boundary  conditions 
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Appendix  B 


Relationship  between  k-f  and 


Two  potential  electrodes 
are  making  contact  with  the 
ground  as  shown  in  the  figure.. 
Prom  the  relation 
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T  =  -  or  vf 

we  concluded  immediately  that  in  the  wire 
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where  subscript  £  refers  to  the  electrode  system. 
But  we,  have 


where  CJU.  is  the  electric  potential  contribution  due 
to  chemical  activity.  We  already  know  that  the  electric 
potential  (fi  is  continuous  across  the  electrode 
contact  (  _  (j£  where  £  refers  to  the  ground), 

therefore  *  ^ 

(fe'-f.-J  =  (ft-fc)  4-  (C.yli-CyU;) 

It  isjobvious  that  Af  and  are  equal  only  when  the 

product  CJU.  is  the  same  for  both  pots.  This  conditions 
is  seldom  satisfied  in  practice  for  the  following  two 
reasons:  (a)  the  lateral  variations  in  the  chemical 
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conditions  of  the  surface  soil  can  produce  drastically 
different  chemical  reactions  at  the  two  electrodes; 

(b)  differences  in  concentration  of  the  non-polarizing 
electrolyte  in  the  two  pots  as  well  as  different 
porosities  can  produce  small  potential  differences  of  • 
a  few  millivolts.  The  second  effect  is  recognized  by 
geophysicists  using  S.  P.  method  and  is  compensated 
for  by  a  "pot  correction." 
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P.  34 

equation  2.4l 
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line  10 
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equation  3.17 
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first  equation 
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last  line 
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line  19 
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"depth  of  d" 

to  "dent.. 

P.  50 

equation  3.40 
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line  21 
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to  ( YJ 

P.  52 

equation  3.24 
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P.  54 

equation  3.30 
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P.  56 
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equation  3.36a 

change 
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equation  5.10.1 
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last  line 
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last  line 
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